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^ , Abstract. We review the properties of the Kronecker (direct, or tensor) product of square 

D I matrices A ^ B ^ C ■ ■ ■ in terms of Hubbard operators. In its simplest form, a Hubbard 

' operator X^^ can be expressed as the n-square matrix which has entry 1 in position 

and zero in aU other entries. The algebra and group properties of the observables that 
define a multipartite quantum system are notably straightforward in such a framework. In 
' (—( ', particular, we use the Kronecker product in Hubbard notation to get the Clebsch-Gordan 

QhI decomposition of the product group SU{2) x SU{2). Finally, the n-dimensional irreducible 

r~| ' representations so obtained are used to derive closed forms of the Clebsch-Gordan coefficients 

' that rule the addition of angular momenta. Our results can be further developed in many 

different directions. 
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1 Introduction 

The Kronecker product, represented by the symbol ®, has attracted the attention of 
researchers in diverse areas of mathematics and theoretical physics over the last decades 
P-[9]. Introduced by Zehfuss in 1858 (see the historical review given in [3]), this is a 
matrix operation also known as direct or tensor product, defined for matrices A = [aij] 
and B of any order to he A ^ B = [aijB]. That is, for instance 



A^B 



Oil ai2 ai3 

0-21 0,22 0,23 



B 



aiiB 


duB 


aisB 


a2iB 


022-8 


023-B 



Matrix calculus includes the derivatives of a matrix with respect to a scalar, a scalar 
with respect to a matrix, and a matrix with respect to a matrix; all these operations 
are defined in terms of the product ® [SHT] (other interesting applications can be found 
in [9]). In physics, this product arises quite naturally if the group properties of the 
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dynamical variables of a given system are considered [T014T2] . and it is fundamental in 
the study of multipartite systems [T3l|ll]. Indeed, the Clebsch-Gordan decomposition 
of the Kronecker product of two irreducible representations is one of the most useful 
problems in group theory, since the reduction of such a product into the sum of irreducible 
representations confirms the unicity of the representation for the simple reducible groups 
[T5] . It is then quite natural to find immediate applications of both, the product ® and 
the Clebsch-Gordan decomposition, in the addition of angular momenta [TB] as well as in 
the identification of symmetries in quantum physics [17J . 

Despite the simplicity of its definition, calculating the tensor product A® B ® C ® ■ ■ ■ 
becomes cumbersome for large matrix sizes and/or for a large number of factors. This 
fact is particularly notable in the design of fast Fourier transform algorithms where the 
factorization of the discrete Fourier matrix F„ is relevant. Namely, for n = 2m with 
m G N, the n-square (Fourier) matrix0: 

Fn = [kj], A,=^(^-^)(*-^), z,j = l,2,...,n, w = e'-'/\ (1) 

can be expressed as the product 

where 

Bn=(^j"' ^ , Ifc = diag (l,l, . ^..,l)^ , = diag{l,w,w'^,. . . ,w''~'^), 

fc-times 

and Ilfc is the k x /c-permutation matrix obtained by grouping the odd columns of the 
identity 1^ first, and the even columns second. Hereafter is the matrix transpose of 
A. The procedure can be repeated if m is even. Indeed, if n = 2*, t G N, then F„ can be 
factorized into the product of t = log2Ti matrix factors (Cooley-Tukey factorization): 

F„ = (Ii ® 5„) {h 8) (I4 8) 5„/4) ■ • • (I„/2 ® 52)Pj, (2) 

with the bit-reversing permutation matrix [18j. Note that each factor 1^ (g) Bn/k in 
(|2]) has only two nonzero entries per row. Thus, only 2n of the entries associated to 
Ifc Bn/k are different from zero. This fact suggests that there must be a better and 
simpler form of calculating the Kronecker product in the factorization algorithms. 

On the other hand, the entries of the first row and column of the Fourier matrix ([T]) 
are all equal to unity while the other entries are either ±1 or ±i. Thus, Fn is a dephased, 
complex Hadamard matrix [19]. The relevance of a dephased matrix is that only its lower 
right {n — l)-square submatrix is necessary in the calculations where it is involved. In 
this subject, it can be shown that the product D^HDc brings any Hadamard matrix H 
into the dephased form for a pair of uniquely determined diagonal unitary matrices 

Fourier matrix Fn — Fn/\/n is unitary. Here we also refer to (HI) as Fourier matrix although it is 
a rescaled version of F„. 
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and Dc [19]. The construction can be simplified since two Hadamard matrices, Hi and 
H2-, are equivalent if there exist diagonal unitary matrices Di and D2, and permutation 
matrices Pi and P2, such that [T9t[20]: 

Hi = D1P1H2P2D2. (3) 

According to the former property, Hi is dephased if Di = D,. and D2 = Dc, no matter 
the form of PiH2P2- In this case dephasing is equivalent to the permutation of rows 
and columns defined by Pi and P2. However, it is not easy to verify whether there 
exist such permutations fT^. The problem becomes even more complicated for large size 
matrices since the permutations grow as A^! Then, it is apparent the necessity of a proper 
framework in which the determination of the above described permutations becomes a 
tractable problem. 

In general, the matrix algebra includes algorithms that are fairly complicated and 
cumbersome for either matrices of large sizes or a large number of matrices to operate 
with. It would therefore be desirable to construct a mathematical framework in which the 
problems like those aforementioned are feasible; no matter the number or the size of the 
matrices involved. To get a suitable approach it is useful to consider the operators X^''^ 
introduced by Hubbard in 1964 [21] (see also [22]). These obey the multiplication rule 

X^'^X^'™ = ^jfcX*'"^, (4) 

and have the properties 

{X''^y = X^'\ J2 ^'''^ = ^' [X^'^ X'='"]± = SjkX''"^ ± 6miX'''\ (5) 

k 

the sublabel in [yl,^]-!- stands for either the commutator (— ) or the anticommutator (+) 
of A and B. 

Hubbard operators provide a way to study groups of particles that interact strongly 
one with each other in such a way that a weak interaction between the groups is also 
allowecj^. For example, they are useful in the description of atoms in which Coulomb 
repulsion prevents double-occupancy of a given orbital [23]. In such cases, the strong 
interactions determine the energy of the groups of particles and can be included in the 
Hamiltonian as linear combinations of the Hubbard operators [2T]i22] . This is the situation 
for strongly correlated electrons [23], whether they are in a cavity [25] or in a two-atom 
molecule [26], as well as for double quantum dots [271|28] among other systems. 

Given an n-dimensional vector space Tin with orthonormal basis {\ipk)}k=i^ the Hub- 
bard operators are written in terms of the outer products of the basis elements: 

X;'^-:=|V'.)(^,-|, t,j = l,2,...,n. (6) 

^To distinguish between the particles of different groups, Hubbard used Xpq [21] as well as Xf [22] to 
denote the operator \i,p){i,q\, with i labeling a given group while \i,p) and \i,q) represent two different 
states of a particle in that group. In contrast, we use a sub-label "n" to denote the matrix order of the 
linear representation of the operator X^'^ , as this is done in equation ©. 
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That is, the operator X^'-^ = \ipi){'4'j\ is a representation of X*'-^ in the space T-in- This 
causes a transition from the state to the state \ipi) of the system that is described by 
the vectors in 7{„. In general, any hnear operator O : T-in "Hn can be represented in 
terms of the Hubbard operators 



(7) 



This property plays a central role in what follows since the operators X^i^ are the cor- 
nerstone of our approach. Using the Hubbard's representation ([7]) one can address the 
algebra of square matrices in compact form, no matter the size or the number of the fac- 
tors. Problems like the determination of the permutation matrices fulfilling ([3]) become 
simpler in this notation. Indeed, if the basis vector \il)k) is the n-tuple that has a unity 
in the fc-th position and all other entries equal to zero, the Hubbard operator ([6]) is in its 
simplest form: 



n 



(8) 



/ " \ 

O(i-i)xo-i) : 0(i_i)x(„-j) 


••• 10 ••• 



0(,i_i)x(j-l) '■ 0(„_i)x(n-i) , 

\ / 

where Qsxt is the null matrix of order s x t. In other words, X^'-' is the n-square matrix for 
which all the entries are zero except the one at the i-ih row and the j-th column, where 
it takes the value 1. Such an array of zeros and a single unit is appropriate to operate 
the Kronecker products ® Bn/k of Eq. ([2]) in plain notation; these products are but the 
linear combination of only 2n Hubbard operators. 

This work attempts to stimulate further progress in the applications of the Hubbard 
operators by introducing a useful manner to calculate the Kronecker product. The paper 
is structured in two main parts. First, in Section [2] some basic definitions and the Hubbard 
notation are introduced. In Section 12] we review some of the most important properties of 
the direct product by using the Hubbard operators as the building-blocks of the Kronecker 
algebra of square matrices. Properties just as the composition of permutations are nicely 
worked in Hubbard notation. Some other properties as the Kronecker powers of operators 
A®^ are explicitly developed for their possible application in quantum control of multi- 
partite systems. In Section H] the useful notation of the direct sum of vector spaces and 
the linear representation of groups are also revisited while the Clebsch-Gordan problem 
is stated in general form. The second part of the paper is devoted to the applications. 
We first review the construction of irreducible representations of the SU{2) Lie group in 
Hubbard notation (Section E]), then general expressions are derived for the representation 
of SU{2) X SU{2) in Hubbard notation (Section EJ. In Section [7] the Clebsch-Gordan 
coefficients of the SU{2) x SU{2) Lie group are expressed in definite form by using the 
Hypergeometric function 3F2. We close the paper with some concluding remarks. An 
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appendix is added to analyze some basic properties of the ceiling and floor functions that 
are required along the paper. 

2 The Hubbard framework 

Matrix calculus was developed for square matrices [29], and finds a lot of applications 
in quantum theory where the observables are represented by Hermitian operators. The 
latter are expressed as square matrices according to the representation defined by the 
mensurable physical quantities and the related eigenvectors, see e.g. [30] • Any observable 
O of a quantum multipartite system S = Si + S2 + ■ ■ ■ is defined on the entire Hilbert 
space H = ® 7^2 ® ■ ■ ■ , and can be expressed as matrix Kronecker products of the 
observables Ok belonging to the subsystems Sk [ISHHj. We shall focus on the properties 
of square matrices in the understanding that they give a suitable linear representation of 
the group of observables defining a quantum system. Yet, most of the results we are going 
to derive can be immediately extended to the case of n x m matrices. However, some 
specific properties of square matrices require caution to be promoted to (or they simply 
can not be applied in) the rectangular case. 

2.1 Notation and basic properties 

Let the ket \x) be an element of the vector space K", with K a field which could be either 
M or C. This will be represented as a single-column matrix containing n numbers Xi G K. 
The latter will be indexed from 1 unless otherwise stated. Thus, 

|x) G K" \x) = {xi,X2,. . . ,Xn)'^, Xfc G K. (9) 

The size n of any n-tuple |a;) will be implied whenever \x) be written as a linear combi- 
nation of the orthonormal basis {|e^)}^^]^ of K", otherwise this will be explicitly stated 
if necessary. Here |e^) means the n-tuple which has a unity in the k-th position and all 
other entries equal to zero. The Hermitian transpose of any ket \x) G K" will be 
represented by the bra {x\, this last is also called the dual of \x), defined as 

\xy := {x\ = {xl,xl,...,xl), (10) 

with x\ = Xfc if K = M, and x\. — Xk for IK — C. The symbol z stands for the complex 
conjugate of 2; G C. Note that the basis vectors are real, i.e. \ef)^ = |e^)'^ = (e^| is 
the n-dimensional row vector having 1 in the k-th position and in all other entries. 
Therefore, the inner product between arbitrary basis elements is non-negative 

{e",\e]) = Skj, fc,jG{l,2,...,n}. (11) 

So that any ket \x) G IK" can be also expressed as the linear combination 

n 

la;) = $^a;fc|efc), 2;, = (e,"|x) G K. (12) 

k=l 
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(13) 



In similar form, 

n 

{x\ = Y,^l{ell xJ = (x|e^)GK. 

k=l 

Hence, the inner product between \x) and \y), both arbitrary vectors in K"', is given by 

n n 



fc,£=l 



k=l 



Since ((x||/))''' = the space is Euchdean (Hermitian) with hnear (sesquihnear) 

metric if K = R (K = C) [I5]. In general we shall write K" = Sp{|e")}"^;^ to denote that 
IK*^ is spanned by the orthonormal set {|e"), jcg), . . . , jej^)}, with n G N. In turn, Sp{|e^)} 
will denote the one-dimensional space spanned by the single basis ket |e^). From the inner 
product ([111) one can identify every (x| with a given mapping of into K. The set of all 
these mappings is spanned by the duals of the basis vectors |e)!) and is included in the set 
of all the functional K" — i- K. We write /C" = Sp{(e"|}"^;^ for such a dual vector space. 
Now, using (fT2|) and (fTOj) . the outer product between \x) and \y) yields the dyad 



\x 



ij=l i,j=l 



n 1 



where the "dyadic" operators 



|en(e"|, ^,j = l,2,...,n 



(15) 



(16) 



are represented by the matrices ([8]). The action of \x){y\ on K" produces the transition 
from the ket \y) to In this sense, the matrix array 



\x){y\ 



\Xn J 



{yl yl 



( x\y\ xiyl 
X2y\ X2yl 



\ Xny\ Xny\ 



xiyi \ 
X2yi 



•^nyli / 



is the linear representation of the transition operator \x){y\ in the vector space 
turn, the action of X*'-' on the basis vector |e)!) gives 

Xll'K) = 5,,\e^), 

so that its action on any ket \x) in reads 

and its matrix elements are easily calculated 



[17) 



In 



(18) 



(19) 



(20) 
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It is then clear that X^^ projects K" into the one-dimensional space Sp{|e")}. The algebra 
of these operators is defined by the inner product, which can be set to coincide with the 
conventional matrix product, but it is simpler to use the algebraic rule 

X'n'X'/ = \e^){e]\e-,){e^\ = 5,,Xli'. (21) 

From ( 1211) . the following result is immediate 

[x;'^ x^-]± = x^'^x^- ± x„^'-x;'^- = 8,,xir ± 5™x„^'^-. (22) 

On the other hand, the operators X^i^ are noninvertible real matrices (i.e., detX^'-' = 0, 

and Xn-' = X^^) such that their transpose and adjoint (conjugate transpose) versions 
coincide 

{Xify = {Xiff = X^\ (23) 
This last expression is easily verified by using ( fT6l) as follows 

{Xiff = {\e7){e]\r = {{e]\n\e'l)r = \e]){e-\ = X^. 

Note that the symmetric operators X^'* are Hermitian and satisfy the completeness rela- 
tion 

n 

I„ = J]X;'\ (24) 

i=l 

Then, the matrices X^i^ correspond to the linear representation of the Hubbard operators 
in the space K", as all the properties ([5]) are verified. 

The introduction of Hubbard operators in the algebra of operators representing quan- 
tum dynamical variables is very advantageous since substantial simplifications are achieved 
with the properties (l2Tl) - (!2^ . Concrete realizations will be presented in the next sections, 
with special emphasis in the square matrix representation. Before that, some words con- 
cerning the case of rectangular matrices are necessary. 

2.2 Rectangular matrices 

To generalize the results of the previous section one would consider rectangular matrices. 
In contrast with the square matrices, a rectangular matrix transforms a vector in the 
space 1-Ln into a vector in the space "Hm where, in general, "H^ ^ 'Hm- It is also well 
known that the multiplication of two rectangular matrices is defined only if the amount of 
columns of the first factor is equal to the number of rows of the second factor. Concerning 
the equivalent of the Hubbard operators in the rectangular case, let -E^xm be the n x m- 
elementary matrix having entry 1 in position (^, j) and all other entries equal to zero [5]. 
This can be expressed as 

Ei;i^ = \e-){e]^\. (25) 
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Note that (Kxm)^ = (Kxm)^ = E^mxn, SO (Kxm)^ and Kxm are defined to act on 
different vector spaces for n ^ m, no matter tlie values of i and j. In tlie same context, 
tlie product between n x m-elementary matrices is constrained to tlie multiplication rule 

rpij rpk,£ _ r pi,£ (0P,\ 

Thus, expressions like E'^^pElf^^ are meaningless since the number of columns of -E^xp 
differs from the number of rows of -E^xm- spite of these apparent complications, the 
theorems of matrix calculus deduced for square matrices may be modified for the rectan- 
gular case. This is particularly right for the "square matrices in the broader sense" defined 
in |Tl], Ch. 2 (see also general expressions in |5j). As we have indicated, our interest is 
addressed to n-square matrices since they represent the most general linear operators in 
the vector space K". The outline above is to stress that care must be taken with regard 
to the generalizations of our results to the rectangular case. 

2.3 Square and permutation matrices 

Using (fT6l) it is easy to see that any n-square matrix A = [aij] can be expressed in terms 
of the Hubbard operators 

n 

A=J2 ('^,J^n', e K, (27) 

so that the conventional matrix product AB is expressed as follows 



n In 



AB=\y^ a,,Xlf ^^^^^n' = E E ^^^^^^^^ = (28) 

\i,i=l / \fc/=l / 1,^=1 V A: / 

where C is the n-square matrix 

n n 

C = E C^/^n'^ = E '^i^khk,t- (29) 

i,l=l k=l 

The complex conjugate A, the transpose , and the adjoint of a matrix A read as 

n n n 

A=Y. a„■X;'^ A^ = ^ a,,X^\ A = ^ 4^^X^\ (30) 

i,j=l i,j=l i,j=l 

On the other hand, the action of A on the basis vectors |e") is derived from fl27|) and f|T8|) . 
this gives 

n 

A\e]) = Y,^kJel). (31) 

k=l 

Then, for an arbitrary vector \x) in K" we have 

n n 

A\x) = J2 (^Kj^^^n^K) = E ^k,eXeK). (32) 

k,j,e=l k,i=l 
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The expression for the trace of a matrix is easily recovered: 

n n n 

(e^Alep = <^kA.kh, = a„- TiA = Y,{e7\A\e1) = Xl^^- (33) 

k/=l i=l i=l 

To give an example, consider a square matrix H = [hij] of size n consisting of unimodular 
entries, \hij\ = 1, and fulfilling HH^ = nl„. This is called a Hadamard matrixlfl. In the 
simplest case, with n = 2 and real entries, we have 

where the factor 1 / a/2 has been introduced to make H unitary. Using (J32l) the action of 
H on \x) & reads as 

H\x) = ^ E = ^ [(^1 + ^2)|e?) + (xi - x,)\el)] . (35) 



From fl28l) . the multiplication of H with itself gives 



1 _^ ^«+^ 



2 

i/=l 



TT 



As another example of interest consider a permutation defined by the bijection tt of the 
set of natural numbers = {1, . . . , n} onto itself. In the Cauchy's two-line notation this 
map reads as 

'1 2 ■■■ n 
7r(l) 7r(2) ■■■ Tx{n) 

In particular, the identical permutation vTe is such that 7re(A;) = /c for all k in S*. The 
set of all n! permutations of S forms the symmetric (or permutation) group Sn of order 
n with the identity TTg as the unit element and the composition of maps as the product. 
Such a group plays an important role in quantum physics (see e.g. Ch. 13 of Ref. [IT] 
and Ref. [H]). For example, the Schrodinger equation is invariant under the permutation 
of electrons since the physical equivalence of all these particles. A linear representation 
of Sn is obtained by assigning a matrix P^r per each permutation vr. This is a square 
matrix of order n that has only one entry 1 per row and column, and is zero elsewhere. 
In Hubbard notation the matrix P^^ reads in simple form 



P. = 5^Xr(^). (37) 



=1 



■^More precisely, the square matrices with ±1 entries and having pairwise orthogonal rows are named 
after Hadamard [3Tj . These are included in the set of self-reciprocal matrices introduced by Sylvester [32| . 
Then, the definition above corresponds to a generalization of what is commonly known as a Hadamard 
matrix |19il33| . 
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In this representation the properties of permutation matrices can be studied in compact 
form. To give a pair of examples consider first the product of and P^r, two permutation 
matrices of order n. From (12T]) we have 

n n n 

k.e=i k.e=i k=i 

Thus, the composition vr o a of permutations vr and a is obtained from the product of 
the corresponding matrices. It is clear that the product of permutation matrices is non- 
commutative as PTrPa = Paow aud aonj^noam general. As a second example let us 
verify the orthogonality of permutation matrices P^^ = Pj^-i = Pj. From (130 p and f l2ip 
one arrives at 

n n 

P^PJ = SnikUi^)Xy = ^n' = (39) 

k,e=l k=l 

Similarly, PjPn = In- From these results it follows the rule (P^Pt^)^^ = P^^P^^. 

To close this section we emphasize that the action of P^^ on any ket \x) G is 
immediately calculated in Hubbard representation 

n n n 

= Y ^^^n^^M) = E = E^-0-)|e")- (40) 

j,k=l j,k=l j=l 

In the next sections some of the properties of the Kronecker product of permutation 
matrices are going to be discussed. 



3 Kronecker algebra in Hubbard representation 

We start the analysis of the Kronecker algebra with the definition of the direct product. 

Definition Kl. Let A = [aij] and B = [fo^.s] be respectively matrices of order 
m X n and k x i over the field K. The Kronecker product A ^ B is the matrix of 
order mk x ni over the field IK defined a.s A^ B = [aijB]. 

As a first example consider the basis vectors |e^), these are matrices of order 1 x n so that 
kfc) ® ) is a matrix of order l-lxn-n = lx rP. Moreover, this n^-tuple has only one 
unity at {k — l)n + j , and is zero in all other entries. In general, the Kronecker product 
of two basis vectors belonging to different spaces |e"^^) and \e^^) gives a tuple of size nin2 
that has a single unit among nin2 — 1 zeros, as this is stated in the following definition. 

Definition K2. Let \e^^) and \e^^) be basis vectors of K"^ and K"^ respectively. 
Then the Kronecker product |e"^^) ® \e"^) is the nin2-tuple having 1 at {ii — I)n2 + i2, 
and zero in all other entries. 
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It is immediate to verify that the vectors constructed according to Definition K2 are 
orthonormal and that there are only nin2 of them. Thus, all of them represent an or- 
thonormal basis of the vector space W^^"-^. We have the next proposition without giving 
a proof. 

Proposition KO. Let = Sp {|e"/)}"^'^^ and K"^ = Sp {|e"^^)}"^'^^ be vector 
spaces. The set of all the Kronecker products |e"^^) ® le"^^) is orthonormal and 
spans a vector space of dimension nin2, written K"!"^ = Sp{|e^^^) ® le"^^), H = 
1, . . . , ni; ^2 = 1, • • • , ^2}, with the following axioms (a, (3, 7, rj are elements of K): 

(i) (c^K)) ® |e-) = a(|e-) ® |e-)) = |e-) ® (a|e-)) 

(ii) («|er;) + (3K)) ® lei;^) = «(|e-) ® |e-)) + /3(|e-) ® |e-)) 

(iii) |er,^) ® ilK) + VK)) = liK) ® K)) + viK) ® K)) 

An arbitrary vector |x) G K"^*^^ can be written either as a twice- indexed linear combina- 
tion 

n=l 12=1 

or as a single-indexed linear combination 

nin2 
fc=l 

Whenever this produces no confusion we shall write le^^*^^"^'"), A; = 1, 2, ... , nin2n^ • • • , to 
represent the basis vectors \e^^) ® \e^^) ® \e^^) ■ ■ ■ , Z£ G {1, . . . , n^}, that span ]K'^i"-2n3-_ 
In particular, if = n for all £ G {1, . . . then 

K"'' = Sp{|er,)®---®|er;r (41) 

is the space of contravariant tensors of rank p while its dual KT'^ is the space of covariant 
tensors of rank p [15] . For instance, in quantum computing it is customary to write |0) 
and |1) for the basis vectors of K^; using our notation they are \e\) and respectively. 
According to Proposition KO, the four products \e\) ® je^^); ^f- ^ {1)2}, span K^. These 
can be expressed in binary form as follows 

\e\) = \el) ® \el) = |0) ® |0) = |00), |et) = \el) ® \el) = |0) ® |1) = |01), 

\4) = \4) ® \4) = |1) ® |0) = |10), \et) = \el) ® \el) = |1) ® |1) = |11). 

The same notation holds for an arbitrary number of two-dimensional factors. For instance, 

|ef) = \el)®---®\el) = |00000). 

^ V ' 

5 times 
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3.1 Kronecker algebra of Hubbard operators 



The Kronecker product A ® B is simple if the factors are Hubbard operators ([8]). In 
this case, the resulting matrix has only one entry different from zero since each of the 
factors has a unique non zero entry. That is, the Kronecker product is closed in the set 
of Hubbard operators. 

Proposition Kl. Let and X^'^ be two Hubbard operators of order n and 
m respectively. The Kronecker product X^ ® X^'^ is the mn-Hubbard operator 
That is, 



m ^ n 



(42) 



Proof. The proof follows from Definition Kl, explicitly 



\ 

O(j-i)xo-i) : 0(j_i)x(„-j) 


••• 10 ••• 





0(„-i)xO-i) 







(n-i)x(n-j) 



( 

0[ra(i-l)+fc-l]x[n(j-l)+f-l] • 0[„(i_i)+fc_i]x[mn-n(j-l)- 



■■■ 10 









\mn-n(i-l)-k\ X [nO-l)+£-ll 







[?Tin-n(i-l)-fc]x[n(j-l)+£-l] 



_ ^n(i-l)+k,n{j-l)+i 

Besides the basic property introduced in Proposition Kl, the Kronecker algebra of the 
Hubbard operators includes the following set of properties. 

Proposition K2. Let X^'"' be Hubbard operators of order a and take A G K. Then 

i) X'^^ ® X^'^ ^ Xy ® X^J' in general. 

ii) {xi;^ ® X'nf = (Xii^f ® {XH'f . 

iii) (XX^^^) ® Xjl'' = X {Xi;^ ® X„^'^) = X^ ® {XX^'') . 
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iv) ixii + x;^^) ® x^J = x^ ® x^/ + x;^^ ® x^-^ 

v) xy ® {x^ + x;^^) = x^/ ® x^ + xl^' ® x:^ 

vi) (x^f ® x„^'0 ® x;'^ = x\^ ® [xl^^ ® x;'^) . 

Proof. Parts iii, iv and v are immediate from Definition Kl. 

i) From Proposition Kl we see that X^-' ® X^'^ = X^'^ ® X^^ requires the roots of the 
system 

n{i-l) + k = m{k - 1) + i, n{j - 1) + ^ = m{£ - 1) + j, 

which in general has no solution for arbitrary fixed values of n and m. A particular 
solution is obtained if n = m, for which one gets i = k and j = i. The symmetric case 
X^ ® X^ is recovered from this last result. 

ii) From (ji2l) and (|23|), 

= x^^ ® xi;^ = {xi;if ® {xyf . 

Remark that this result immediately gives (X^^ ® X^'^)^ = (X^^y ® (X^'^)^, since the 
Hubbard operators are real matrices. 

vi) From fl42D . 

{xi;i ® xy) ® x;-^ 

_ ynii-l)+k,n{j-l)+e yr,s 

— ^mn yy 

_ -^p[n(i-l)+fc-l]+r,p[n(i-l)+f-l]+s _ ^pn{i~l)+p(k-l)+r,pn{j-l)+p{i-l)+s 

— -^mnp — -^mnp 

= x\^ ® X? = x^ ® {xy ® x;'') . □ 

Proposition K2 includes the basic properties of the Kronecker product of Hubbard opera- 
tors. Applied to the set {Xf '"^j, they mean that the product ® is distributive over ordinary 
matrix addition (iv, v), associative (vi), compatible with ordinary matrix transposition 
(ii) as well as with matrix multiplication by an scalar (iii) and, in general, non-abelian (i). 
Next we generalize such properties to the case of arbitrary square matrices while some 
other algebraic relationships are derived. 

3.2 Kronecker algebra of permutation matrices 

The Kronecker algebra of the Hubbard operators is particularly useful in the operating 
with, and the construction of permutation matrices. In this section we report some 
results on permutation matrices that are fundamental in the ensuing applications of the 
Kronecker product. From Proposition Kl one has the following result. 
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Theorem PI. The square matrix 



n = ^jQ^^xr 

is a permutation matrix of order n^, defined by the rule 

7r(p) = n{p + n — 1) — {v? — p = 1, 2, . . . , n^, 

with 



(43) 



P 



(44) 
(45) 



the ceihng function apphed on ^ (see Eq. (]A.2|) of the appendix). 
Proof. From fj42|) we have 



Let us define p = n{i — 1) + j, then p = 1,2, ... ,ni, and j = p — n{i — 1). Therefore, 



n m 



n = E E K 

i=l p=n(j-l)+l 



p, n(p+n— 1) — (n^ — l)i 
2 ) 



with 



p » In 

-<i<- + l- -<- + l. 
n n n n 



(46) 



Using Lemma Al(i) of the appendix we get 



to write 



^p, n{p+n—l) — {n^ — l)p' 
o2 



p=i 

Comparing this last result with (l37|l we arrive at the rule (144|) . It can be verified that 
tt{p) is indeed a bijection on the set {1, . . . , n^}. □ 



To get some insight on the meaning of the permutation matrix (1431) let us consider an 
arbitrary contravariant tensor of rank 2: 



\Xi 



(47) 



*1,«2 = 1 



il,*2 = l 
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The action of 11 on this last vector gives 



\i,j=^ / \i, 1*2 = 1 

n 



i:i,*li«2 = l «1,«2 = 1 

Thus, relative to the indices labeling the contravariant tensor fH7|) , the operator 11 corre- 
sponds to the bijection 112 : (1,2) 1— )■ (2, 1). Hence 11 = G 5*2. Indeed, there are only 
2! = 2 different permutations on the set {1, 2}, these are the identity vTe = tti and 7r2. In 
Hubbard representation we have 

n n 

P., = 5^ ® Xi;^ ^ Yl ^^^^n ® P^. = n. (48) 

i=l *ii=l 

From these results it is easy to verify that 

5(2) = |(P., + P.J = ^ E (1 + ^*.)^n' ® ^ii' (49) 



is the symmetrization operator for the vectors in K" . Namely, 

C !\ \^\ \\ ® \X2) + \X2) ® \Xi) 
S(2){\Xl) ® \X2)) = 

is a symmetric tensor of rank 2. In a similar form, the operator 



1 " 

A2) = I [x(vri)P., + X{^2)P.,] = 2 E [^(^1) + ^(^2)5.,] X^'^- ® (50) 
with x(^) the parity of the bijection vr [15j, produces antisymmetric tensors of rank 2: 

, I \Xi) ® \X2) - \X2) ® \Xi) 

The generalization of the above results to tensors of arbitrary rank is straightforward. We 
summarize this in the following proposition without a proof. 

Proposition PI. The operators 

1 ^ 1 ^ 

^' e=i ^' i=i 

with TTg G Sp, TTi = TTe, and P^^ a definite linear combination of the Kronecker 
products 

Xn,n ^ xi^'^^ ® ■ ■ ■ X;-^^ Ik, jfc e {1, . . . , n}, 

produce respectively the symmetrization and antisymmetrization of the contravari- 
ant tensors of rank p. 
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The operators S^p) and A{p) are useful in group theory and symmetries [T0lfTT| [T5 t [T7]. 
In the hterature of combinatorics they appear in connection with the concepts of deter- 
minant and permanent of a matrix, these last give rise to entire treatises ^34j and are 
associated to the concept of majorization of vectors that is fundamental in the algorithms 
of quantum computing [H] and in the geometry properties of the quantum states [35] 
as well. The symmetrization of vectors as this has been indicated above is also useful in 
the analysis of the Majorana representation of multi-qubit states for the studies of the 
barycentric measure of quantum entanglement [36|- The major result in Proposition PI is 
that the Kronecker products defined in Theorem PI as permutation matrices are noting 
but the building blocks of the symmetrization and antisimmetrization operators of the 
contravariant vector space fHTj) . The result reported in Theorem PI has been already 
included in the works by other authors (see e.g. Eq. (4) of Ref. [T], and Section 2.5 of 
Ref. [3]) but, as far as we know, such works give no reference to the explicit form of the 
permutation. Here, equation ( H^ gives the concrete realization of such a permutation 
and the connection of 11 with the operators (150|) has been also achieved. 

On the other hand, the Kronecker product of permutation matrices is also compatible 
with the composition of permutations described by equation That is, the product 

® is closed in the set of permutation matrices. 

Theorem P2. Let Pn{n) and Pa{m) be the n and m-permutation matrices defined 
by the rules vr and a respectively. The Kronecker product P-jrin) ® Pai"^) is the 
nm-permutation matrix Pa{n,m) defined by the rule 

a{p) = m[7r(p') — 1] + cr(p — mp' + m), (52) 

with p' = \^] . 

Proof. From Proposition Kl and the linearity of ® one gets 

n m 

Pa{n,m) = P^{n) ® P,(m) = ^ ^ 

i=i j=i 

Following the proof of Theorem PI we realize that the change p = m{i — 1) + j gives 
rise to equation fHBl) with ^ m, so that p' = \^~\ and the rule ( l52i) follows from the 
expression 

nm 

p=i 

Now, using the multiplication rule ( 12T1) and property ( l30l) we have 

nm nm 

P^{n,m)P^{n,m) = ^ Xl^)'^^^^'^) = Y.^n^^'"^'^ = ^nm, 

p,q=l p=l 

and a similar procedure shows that Pa{n,m)P^{n,m) = I^m- ^ 
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Now, let us consider Proposition K2(i). This indicates that the Kronecker product of two 
Hubbard operators, X^i^ and X'^^, is non-abehan in general. Such restriction, however, can 
be relaxed because Xli^®X^^ has only one entry 1, as this is established in Proposition Kl; 
the same is true for X^^ ® ^n^- Therefore, it should be possible to arrive at X^'^ ® XJ^^ 
by applying the appropriate permutation of rows and columns in X^i^ ® ^m^- That is, 
X^^ (g) X^^ and (g) X^'-' must be permutation equivalent. 

Proposition P2. The Kronecker product X^'-' ®X^'^ is permutation equivalent to 
X^^ ® X^^ . That is, there exist P^r, a permutation matrix of order nm, such that 

Pj {Xi;^ ® X^') P. = X^^ ® X;'^-. (53) 

Proof. We use the multiplication rule f l?I]) and the expression together with the 

linearity of ®, to arrive at 

V*>i (<7i V'^i^^ P yir(mi—m+k),n(mj—m+e) 

Hence, in order to satisfy fl5^ we have 

7r(mz — m + k) = nk — n + i. (54) 

The bijection vr we are looking for is defined by this last equation (the labels j and i 
satisfy the same equation under the change i j and k i). □ 

The permutation equivalence of matrix Kronecker products is of enormous interest 
in quantum information theory as this is useful in solving the problem of constructing 
maximally entangled bases of multipartite quantum systems. If dealing with Fourier 
matrices, it is possible to discriminate whether the permutation equivalence is preserved 
for the Kronecker products if matrix multiplication by unitary diagonal matrices is also 
allowed [19]. As discussed in the introduction, this problem deals with equation (|3]) for 
which the permutation matrices Pi and P2 are to be determined. The results we have 
presented in this section are addressed to show the basic operation rules of the permutation 
matrices when they are expressed as linear combinations of Hubbard operators. Further 
insights will be given in the sequel. 

3.3 Kronecker algebra of square matrices 

We now consider the properties of n-square matrices associated with the Kronecker prod- 
uct. One way to phrase the main subject of this section is to say that every n-square 
matrix A is a linear combination of Hubbard operators, just as this is stated in Eq. (l27|l . In 
such a representation the assertions of the following propositions, theorems and corollaries 
are readily verified. 
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Theorem Ml. The Kronecker product of A = [aij] and B = [bk/], respectively n 
and m-square matrices, can be written as 

nm 

That is, the Kronecker product A ^ B is a hnear combination of the Hubbard 
operators of order nm. 

Proof. From Proposition Kl and the linearity of we have 

n m 

i,j=l k,e=i 

Let us define p = m{i — 1) + k and q = m{j — 1) + i, so that p,q = 1, . . . mn. Therefore 
k = p — m{i — 1) and £ = q — m{j — 1). Hence, 

n ni mj 

A (8) B — ^ ^ ^ ] ^ ^ (^i,ji'p^rn-mi,g+m-mj-^mn- 

i,j=l p=m{i—l)+l q=m{j—l)+l 

The indexes i,p and m in the third sum are related by Eq. ( H6|) . The same is true for 
the indexes j,q and m in the fourth sum, so that according with Lemma Al(i) of the 
appendix we get 

^"^ q'. (55) 



p_ 
m 



P, j 



m 



Then 

nm 

A (S) B = ^ ^ Cp^qX^^, Cp^q '.= ttp'^q' hpj^rn—mp' , q+m—mq' ■ I— ' (^6) 

p,q=l 

As a first consequence of this theorem we realize that A = I„ and B = l^n produce 
C = Inm- That is, the Kronecker product of identity matrices is an identity matrix. 

Corollary Ml.l. The identity is preserved in the Kronecker product. That is 
Proof. Use Theorem Ml with aij = 6ij and b^/ = Ske. □ 

Another important consequence of Theorem Ml is that the nontrivial Kronecker powers 
of matrix A, written A'^^'^^ with k > 1, also admit a definite expression in Hubbard 
notation. 

Corollary Ml. 2. The Kronecker product of the n-square matrix A = [aij] with 
itself k > 1 times, denoted is given by the expression 

A^'+'=J2<^iT^K'.., k>h (57) 

p,q=l 
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with 



and 



p,<i 



fc-i 

n 



^Pk,gk XI ^Ps+n-nps+i,qs+n-nqs+iy 
s=0 



Ps 



p_ 



s = 0,1, . . . ,k. 



(58) 
(59) 



Proof. From Theorem Ml, with B = A in fl56|) . one has 

^2 



^ ^ ^ (^p' ,q' (^p+n—np' , q+n—nq' -^^2 ; 



p,i}=l 



with and g' as they have been introduced in ( l55l) . Applying again Theorem Ml with 

A o A'^^ and B ^ Awe get 



.„3 



® A — [a(^piy^(^gyapi+n-n{p'y,q'+n-n{q'y]o.p+n-np',q+n-nq'X'^2 (60) 



where 



p,ij=i 



'x'' 










n 




n 







, x = p,q. 



In the last result we have used Lemma Al(iv) of the appendix. Here, it is convenient to 
use the notation introduced in (1591) . so that (!60l) reads as follows 



A^ 







^ ^ ^P2,q2 




p,q=l 








^ ^ ^P2,q2 


(n 


p,q=l 


\s=0 







p-.q 



p,q=l 



P,<1 
2 ■ 



The proof is completed by induction. □ 



The action of A*^* on the vector space Sp {|e")®*}"^-,^ represents the parallel action of 
t operators A on t vector states \ip) G Sp{|e")}"^^. This manifestation of the quantum 
parallelism is a fundamental feature of many quantum algorithms [H]. Thus, quantum 
circuits can be constructed to evaluate a function f{x) for multiple values of x simultane- 
ously. Most of the procedures implemented to calculate functions on an arbitrary number 
of bits use the Hadamard transform if®". This operation is just n Hadamard operators 
acting in parallel on n qubits. It is then profitable to get a practical expression of H'^" 
particular application of Corollary Ml. 2. 
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Proposition Ml.l. Let H be the Hadamard matrix ( IMIl . then 



2fe+i 



p,g=l 



with Ps = \^] , = , and 



P- Q ■■= ^{Ps - - !)• 



s=0 



Proof. For n = 2 and A = H, equations ( 157|) and (!58|) respectively read as 

2fc+i 

and 



fc-1 



7,(fe+i) _ /, TT /, 

'Ve "'Pk,qk \_\_"'Ps-2-2ps+i,qs-'2-2qs+ii 

s=0 



X 

2^ 



X = p,q. 



From (IMl) we know that 



Therefore, 



''Ps-2-2ps+i,qs-2-2qs+i 



(^_l)(P«-2-2ps+i-l)('2s-2-2g^,+i-l) ]_~) {Pa-l)fe-l) 



and 



/^(fc+i) = ^ 



4)ELo(p--i)('?=-i). 



(61) 



(62) 



(63) 



Equation fl6T]) follows from the introduction of this last result in f l63|) . □ 

Let if„ denote a Hadamard matrix of order n; the matrix (= H2) defined in Eq. (!34l) 
and used in (|6Tll is the simplest example. The next case is found for n = 4 as the 
orthogonality condition on the rows of if„ forces n to be even |33]. Proposition Ml.l. 
gives an easy way to construct a Hadamard matrix of any order because H®^^^ is of size 
2fc+i_ Yoi instance, H4 = H ® H = if®^ reads as follows 



1 ^ 1 

= i7®2 = - ^ (_i)(p-i){g-i)+(rii-i)(r§i-i)xP''' = - 



p,g=l 



/I 1 1 1 \ 

1-1 1-1 

1 1-1-1 

VI -1 -1 1/ 
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The matrices if„ are known as symmetric multiports (or Zeilinger matrices) in quan- 
tum optics and have apphcations in combinatorial problems, coding algorithms and 
quantum engineering, among a diversity of subjects. 

In order to appreciate the significance of Proposition Ml.l., let us apply the operator 
ff(^k+i Qj-^ ^^-y q£ ^Yie vectors spanning . Using ( 16T|) and (I3T]) with A = H®'^^^ we get 

^«^+i|ef-) = 1 '£(_i)Eto(..-i)0.-i)|ef-), (64) 

with Us = [^1 for y = p,j. In particular, if /c = 1 the latter expression gives 

^®2|g4^ = 1 ^(_i)lP-i)U-iHlPi-i)ai-i)|e4^. (65) 



Explicitly, 





(p-i)(j- 




p=l 






' |ef) + 


|e^) + 


|e^) + 1 




|e^) + 


|e^) - 1 


|ef) + 




14) - 1 


. |ef) - 




|e^) + 1 



,4\ 



el), J =2 
64), J =3 
el), J =4 

We can see that the action of H®'^ on |ef) produces an equal superposition of all basis 
states. In quantum computing, this corresponds to set an empty quantum register of 
2 qubits |00) = |ef) into an equally weighted distribution of all the basis states of the 
register |00), |01) = |10) = jcg) and |11) = le^). At this stage, it would be useful to 
show the translation of the results from our notation to the binary one, which is widely 
used in the quantum computing context. We first give the following definition. 

Definition Ml.l. Consider a positive integer x < 2'=+! with keN. The expansion 
of X in powers of 2 is defined by the binary coefficients Xs G {0, 1}, s = 0, 1, . . . , fc + l, 
as follows 



k+l 

X = 

i=0 



^x,2\ (66) 



Now we pay attention to the coefficients of the linear combination ( l6Ti) . as only these 
must be rewritten to get equations ( 16T|) and ( IMl) in binary form. The next proposition is 
necessary. 

Proposition Ml. 2. Let p and q be respectively the i-th and j-th powers of 2 with 
i, j = 0, 1, . . . , A; + 1, and A; e N. Then 

(_l)Eto(r*l-i)(r*l-i) = (_i)ELo(p-i)49-i)=^ (67) 

where {p — l)s and (g — 1)^ are the s-th binary coefficients of p — 1 and q — 1 
respectively. 
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Proof. We first prove the identity 



(68) 



with \_x\ the floor function of x (see Eq. ( ]A.3I) of the appendix). Using the binary 
expansion of p and q we can write 





P 


-1 








Given s one has 


[x,2' 





fc+l 



(69) 



i,j=0 



•,e-s 



p 






.2'. 




.2-'. 



all the terms labelled with either i < s or j < s in (169|1 are equal to zero. We have four 
partial sums 

fc+i fc+i fc+i 

j=s+l i=s+l i,j=s+l 

The second and third terms of this last result are either zero or an even number, so they 
can be omitted from the exponent of —1 in the expression at the left of Eq. (l68l) . The 
same can be said of the fourth term after the change i — s — )■ r, j — s — > t. Finally, 
[PsJ IqsI = Psls since Psils ^ {0, 1}. Therefore, taking into account only the elements of 
that contribute in nontrivial form to the exponent of —1 we get (1681) . That is. 



p 




Q 


* 


.2'. 




.2'. 





Pf.qg ^ ( — l)E3=oL2^J L^J = ^ j^^Es^oP^'Js^ 

Using Lemma Al (part v) of the appendix and this last result, equation (167|) follows. □ 

To verify the compatibility of our results with those obtained in a binary representation 
let us rewrite the coefficients of the linear combination (IMl) according to Proposition Ml. 2. 
For k = 1 one has 



1 4 

H^'^\e^) = - ^(-i)(p-i)o(i-i)o 



+{p-i)iO-i)iU4 



(70) 



which explicitly gives the following result 



' (-1. 


0.0+0.0|g4^ + ( 


_-|^^l-0+0-0 


ei) + 


^_l)0-0+0-l 


el) + { 


-r 


1-0- 




j 


= 1 


(-1 


)0.1+0.0|g4^ + ( 


_iy-i+o-o 


4) + 


'_l)0-l+l-0 


et) + i 


-1 




''■'let) 


j 


= 2 


< 

(-1 


r+°->t) + ( 




et) + 


^_X)0-o+i-i 


et) + { 


-1 


iO-l- 


''■'let) 


J 


= 3 


. (-1 




_]^^M+0-l 


et) + 


._1)0.1+1.1 


et) + { 


-1 


)^-^- 


"■'let) 


j 


= 4 
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The comparison of (TfOjl with ( IMI) shows that the results are consistent in both repre- 
sentations. The final step is to express \e^) in binary form. Using Proposition KO and 
Definition Ml.l we make 



:= |(j - l)o, (j - l)i, . . . , (j - 



with (j — l)s G {0, 1} the binary coefficients of j — 1 up to 2'^. If A; = 1 then \ep — >■ 
\j — = |(j — l)o, (j — Hence \ef) — |0)(i) = |0,0), so that we can write 

\ef) = 1 00), and so on. Then, we can write equation ( l64l) in the standard form 



where n = 2'''^^. Here \x) and \z) are in binary notation. 

Coming back to Theorem Ml, we stress that this can be generalized to an arbitrary 
number of factors by including the products A ® B and A^^^^ as particular cases. This 
is stated in the following proposition. 



Proposition Ml. 3. Let Ar = a^fj be a square matrix of ordern^. The Kronecker 

product y4i (g) (8> ■ ■ ■ ® Afc+i is the square matrix of order n^'^^ = nin2 ■ ■ -n^j^i, 
expressed as the following linear combination of Hubbard operators 



where 



A = Ai® A2 



V,1 



Pk 



k-1 



At+i = ^ a,;," x'j„ k > 1, 

p,q=l 



, {k-s+1) 

'Ps+"fe-s+l-"fe-s+lPs + l. qs+rik-a + l-nk-a + lls + l ' 



(71) 



(72) 



and 



Ps 



p 



(73) 



Proof. This is immediate by following the proof of Theorem Ml and Corollary Ml. 2. □ 

The advantage of having an expression for the matrix Kronecker product as general 
as the one reported in Proposition Ml. 3 relies on the fact that this includes an arbitrary 
number of factors, the size of which is in turn arbitrary. Immediate results can be obtained 

the product (17T!) 



(1) 



(2) 



as particular cases. For instance, if A; = 1 and a 

gives the expression of A® B reported in Theorem Ml. Now, let all factors in (17T|) be of 
the same order, namely = n for r = 1, 2, . . . , /c + 1, then the coefficient (172|) reads as 



a = a 



fe-i 

"'Ps+n-nps+i, qs+n-nqs+i^ 

s=0 



(fc-s+1) 



(74) 
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and the definitions of ps and in ( 1731) are reduced to tfie ones given in ( l59l) . Using these 
last results and considering the situation in which all the factors are equal, i.e. Ar = A 
for all r, one recovers the expression for reported in Corollary Ml. 2. In such a case, 

the super-index could be omitted from all the matrix elements appearing in (17^ . On the 
other hand, let A^ be Fourier matrices Fr, then (I7T1) is the 'Fourier Kronecker product' 
defined in [3H], expanded in terms of Hubbard operators. In addition, this is 'factored' if 
the size of the matrices A^ are natural powers of prime numbers, and is called 'pure 
factored' if such powers are ordered: 

F = F^i, (g) ® ■ ■ ■ ® F^i,+, , ii»i2» ■■■ » 4+1, e N, 

where a is a prime number. The immediate generalization considers rir = a^'' , with 
standing for a prime number and r = 1, . . . , + 1. In Ref. a 'multi-index' notation is 
also introduced to operate the above mentioned products. Such notation works well for 
the necessities indicated by its author in the study of Fourier matrices. However, it could 
be not easy to generalize this to the Kronecker algebra of square matrices other than the 
Fourier ones. In contradistinction, the algebraic rules presented here are based on the 
linear superposition of Hubbard operators, so that they minimize the number of indices 
that are required to operate with. In this form, the Hubbard representation facilitates 
the calculation of the Kronecker algebra of the square matrices of any sort. 

Theorem M2. Let A = [aij],C = [cp^g], and B = [bk/],D = [dr^s], be pairs of 
n and m-square matrices respectively. The usual matrix product of the nm-square 
matrices A (g) 5 and C ® D fulfills 

{A®B){C®D) = AC®BD. (75) 

That is, the Kronecker product of square matrices is compatible with ordinary 
matrix multiplication. 

Proof. From Proposition Kl and the linearity of ® one arrives at 

n m 

i,j k,l 
n m 

p,q r,s 
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Using ( l2Til and Lemma A2 of the appendix we get 



n m 



i,j,P,g k,i,r,s 



i,j,q k,£,s 
n m 



i,q k,s 

' n \ / m \ 

Y}^AC\,X^ ® Y,{BD\sX'^^^ \=AC®BD. □ 

i,q / \ k,s / 

At this stage we have to stress that Theorem M2 can be extended to rectangular matrices 
whenever the involved matrix products make sense (see the discussion of Section 12. 2p . 
This fact will be taken into account in e.g., Corollary TM2.2 and Proposition G.l. 

Corollary TM2.1. Let A and D be square matrices of order n and m respectively. 
Then the following relationship holds 

A®D= {A0lm)iIn®D). (76) 



Proof. Use B = 1^ and C = I„ in Theorem M2. □ 

Corollary TM2.2. Consider the eigenvalue equations A\ai) = ai\ai) and B\bj) = 
/3j\bj), with i = 1,2, ... ,n and j = 1,2, ... ,m. Then (i) the nm numbers OiPj are 
the eigenvalues of A(^ B associated to the vectors \ai) ® \bj) (ii) The eigenvalues of 
A ® Im + In ® -B are the numbers a, + /3j. 

Proof. Using f[T^ the vectors |aj) and \bj) can be represented by row tuples of size n and 
m respectively. Then, Theorem M2 applies as follows 

{A ® B){\ai) ® \bj)) = A\ai) ® B\bj) = ail3j (la^) ® \bj)) , 

and the proof of (i) is completed. In similar form, 

{A(^Im + In ® B){\ai) ® \bj)) = A\a,) ® + I^ja^) ® B\bj) = (a^ + Pj){\a,) ® \bj)) 

completes the proof of (ii). □ 

The following properties are presented here in connection with the Kronecker algebra of 
Hubbard operators, because any square matrix can be expressed as a linear combination 
of such operators. In this form, the proof of each item is simple by using equation (1271) 
and Proposition K2. 



25 



Theorem M3. Let A, B and C be n-square matrices and A G K. Then 

i) In general, A® B ^ B ® A. 

ii) {A ® B)^ = A^ (^B^ and {A ® Bf = ^ B^ . 
hi) (XA) ® B = \{A® B) = A(g) (XB). 

iv) {A + B)®C = A®C + B®C. 

v) A0{B + C) = A®B + A®C. 

vi) {A® B)®C = A®{B ®C). 

Proof. The theorem foUows from Proposition K2 and properties (l27j) and (l30l) . For 
instance, the proof of both parts in (ii) is quite similar: 




1 1 



k,e=i 



Ki,j,k,e=l 




i,j,k,£=l i,j,k,l=l 
= At® fit. 

Then (A ® 5)^ = (AWbY = A^0B'^ = A'^^ B^ . □ 

As for Proposition K2, the above properties mean that the Kronecker product is dis- 
tributive over ordinary matrix addition (iv, v), associative (vi), compatible with ordinary 
matrix transposition (ii) and the matrix multiplication by an scalar (iii). Remarkably, 
although the product ® is non-abelian for arbitrary square matrices (i), it is possible to 
set equivalence classes between Kronecker products that are 'abelian' up to a permutation 
matrix. Our claim is based on the Proposition P2 as well as the linearity of the Hubbard 
operators and leads to the following theorem. 

Theorem M4. Let A = [aij] and B = [bk/] be two square matrices of order n and 
m respectively. The Kronecker product A^ B is permutation equivalent to B ® A. 

Proof. From Proposition P2 we know that there exists a permutation matrix P such that 
( l53l) is true. Then, by linearity in the conventional matrix product we have 

n m 

P^iA ^B)P= EE "^.^•^'^.^ [P"^ i^n ® X'J) P] 

i,j k,£ 
n m 

= EE^M«.,.^:^'®^;''' = i?®A □ 

i,j k,t 
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The next properties show that the Kronecker product is compatible with the conventional 
measure properties of square matrices. 

Proposition Ml. 4. Let A = [oj .,■] and B = [bk/] be two square matrices of order 
n and m respectively. Then 

Tt{A (S)B) = Tr(A)Tr(5) = Tr(5 (g) A). 
Proof. Using Proposition Kl, the linearity of (S), and property (!20|) we can write 

n m 

(e""!^ (g) -Ble"") = ^ ^ aijbk,e 5s,m(i-i)+fe Sm{j-i)+£,s- 

i,j k,£ 

Then, ( 133|) gives 

nm n m 

n m / " 

= 0'i,jbk/ ^ij^k/ = I 

= Tr(A)Tr(5) = Tr(5)Tr(A) = Tr(5 ® A), 
where we have used Lemma A2 of the appendix. □ 

Proposition Ml. 5. Let A = [aij] and B = [bk/] be two n-square matrices. Then 

Det(A O 5) = (Det v4)"(Det 5)". (77) 

Proof. From Corollary TM2.1 we have 

A^B = 

Then 

Det{A(g}B) = Det(A(g)I„)Det(I„(g)5) = Det(A O I„)(Det 5)". 

The last term in the previous expression is because I„ (g) 5 is a Jordan matrix having n 
repetitions of B along the diagonal. Now we use Theorem M4 to get 

Det [P^iA In)P] = Bet{A » I„) = Det(I„ A) = (Det A)". 

Equation fl77p follows from the above results. □ 

Additional properties (and proofs) of the Kronecker product of matrices can be found 
in the books [5HZ]. The most recent summary of the properties of the ® operation has 
been reported in [9] (see also [8]). Although these references are addressed to the case 
of rectangular matrices, most of the useful applications of the product ® require square 
matrices only. Indeed, properties like the permutation equivalence of Kronecker products 
are fundamental in the simplification of algorithms to process data or in the study of 
symmetries associated to a given physical system, among other applications of square 
matrices. Next, we shall apply our results in the analysis of the angular momentum of 
quantum systems. 
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4 Basics of the Clebsch-Gordan decomposition 



In this section we consider the Kronecker product of two irreducible representations of 
a given group. Our interest here is to analyze the generalities of the reduction of such 
a product as a sum of irreducible representations. To start with, wc require some basic 
definitions of the direct sum of vector spaces and the linear representation of groups. 



4.1 Direct sum of vector spaces 

Definition SI. Let 1-1 = 1-1'® H" be the direct sum vector space of H' and H" 
over the field K. Then H' CiH" = Sp{|0)}, with |0) the null vector of H. Assuming 
Dim('H) < oo, we have Dim('H) = Dim('H') + Dim('H"). In general, any vector 
\x) & H can be written in the form of a two- vector tuple \x) = {\x') , \x"))'^ , with 
\x') e H' and \x") e H". The following matrix notation will often be used 



\X) 



/ ^1 \ 

Xn' 
\ ^n'+n" / 



( < \ 



\X 



\X 




X^i 



Now, let A' -.li' ^ %' and A" : %" H" be respectively automorphisms of H' and H". 
Assume their action is as follows 

A'\x') = \y'), A"\x") = \y"). 

They together can be expressed in a single matrix operator acting on the entire space 
H' © H" in block-diagonal form 



A' 


On'xn" 


On"xn' 


A" 





We shall write in this case A = A' (B A", and we shall say that the automorphism A : 
H ^ H decomposes into the direct sum of A' and A". 

Definition S2. Let H = H' ® H" be the direct sum of vector spaces %' and H" ^ 
with Dim('H') = n' and Dim('H") = n" . The direct sum operator A = A' ® A", 
written in diagonal matrix form as 



A = A'®A" = 



A' 


On'xn" 


On"xn' 


A" 


as A\x) 


= {A'e 



with A' e Aut(?^') and A" e Aut(?^") respectively 
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Hereafter Aut(X) denotes the set of all automorphisms of X. The next proposition allows 
to distinguish between operators that can be decomposed as a direct sum when they act 
on any vector space 0^ Tik, and those that act on 0^ Tik in a more complicated manner. 

Proposition SI. The most general operator acting on the direct sum space T-L = 
W © H" is of the form 

A -- 



A' 


B 


C 


A" 



with B and C matrices of order n' x n" and n" x n' which respectively maps %" 
into %' and %' into %" . A' and A!' are n' and n"-square matrices respectively. 



Proof. The proof is immediate in matrix notation 

A\x) = 



A' B \ f \x') \ _ f A'\x') + B\x") 
C A" J \ \x") J ~ \ C\x')+A"\x") 

where the products B\x") and C\x') are n' x 1 and n" x 1 matrices respectively. □ 



4.2 Linear representation of groups 

Definition Gl. Let G and H he a, group and a vector space respectively. The 
homomorphism 

T: G^Aut{H) 
9 ^ T{g) 

is a linear representation of G on "H. The vector space "H is called the representation 
space induced by T. The operator T{g) fulfills 

T{g29i) = T{g2)T{gi), V^i, g2 G G. 

Henceforth all the operators will be expressed in matrix form, so that the representations 
considered will be linear. 

Definition G2. Given a representation T of G on "H, a subspace Tig of "H is said 
to be G-invariant if, for all |a;) G Tis and all g E G, we have T{g)\x) G Us, i.e., 
T{g) G Aut('Hs) for all g & G. In any representation there exist two trivial invariant 
subspaces; H itself and the null vector space Sp{|0)}. A representation T of G on 
v. is irreducible, written A(G), if the only G-invariant subspaces of H are Sp{|0)} 
and H itself. 

Proposition Gl. Let A'(G) and A"(G) be two irreducible representations of 
a given group G. Let Dim('H') = n' and Dim('H") = n" be the dimensions of 
the corresponding representation spaces. Given g G G, the Kronecker product of 
matrices A'(gf) and A"{g) is a square matrix of order n'n" acting on Ti' ® "H". That 
is, the set of matrices 

{A'^A"){g)^A'{g)^A"{g) 

defines a hnear representation of the product G xG, with %' ® H" as the represen- 
tation space. 
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Proof. Given \x') E W and \x") E n", one has A'(^)|x') = \y') E W and M'{g)\x") = 
\y") E H" . Using Theorem M2 we obtain 

{A\g) ® A"ig)) ® \x")) = A'{g)\x') ® A"{g)\x") = \y') ® \y") = \y). 

Therefore, given g E G, there exists a matrix T{g) = A\g) ® A"{g) E Aut('H' ® l-i") such 
that T{g)\x) = \y), with \x) = \x') ® \x") and \y) = \y') (g) \y"). □ 

4.3 Clebsch-Gordan decomposition 

Let us stress that in general the Kronecker product A' (8> A" is not irreducible, even 
though A' and A" are irreducible. This unpleasant situation defines the problem of 
reducing A' ® A" as much as possible to a representation where the vector space H is 
G-invariant. The best we can do is to split the representation space H into a set of 
subspaces Tik fl Tij = Sp{|0)}, k ^ j, each of them being G-invariant. Then, A' (g) A" is 
decomposed into a set of irreducible representations A^, one per each subspace Tik- 

Definition G3. Let A'{G) and A"{G) be two irreducible representations of a given 
group G. The Clebsch-Gordan decomposition of (A'® A") ((7) is the reduction of the 
Kronecker product A'{g) (g) A"{g) into a direct sum of p irreducible representations 
Ak{g) defined as 

p 

{A! ® A!'){g) = Ak{g) = A,{g) © A^{g) © ■ ■ ■ © Ap{g), \fg E G. 

k=l 

Note that the representation space of the Clebsch-Gordan decomposition A^ is the 

direct sum 0^^^^ Tik- Therefore, Definition G3 implies that the vector space Ti' © Ti" 
is decomposed into the direct sum H' © Ti" = ^l^^Hk- If it is possible to solve the 
Clebsch-Gordan decomposition problem for a given product A' © A" we say that such 
representation is completely reducible. Remark that only finite dimensional representa- 
tions are always completely reducible, for infinite dimensional representations this is not 
generally true [TT]. 

5 The group SU{2) in Hubbard representation 

Let us consider the group of unimodular (i.e., with determinant -|-1) unitary 2x2 matrices 
SU{2). The general form of one of these matrices is 

A=(J^). |a|^ + W^ = l. (78) 

In terms of the identity I2 and the Pauli matrices a = {(Tx, Cy, a^), with 
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equation ( 1781) reads 



A = ao^2 + ■ a , a = (oi, 02, 03), a = ao + ia3, b = a2 + iai. 

The well known correspondence between the elements of SU (2) and the points of a sphere 
of radius vr in the three-dimensional euclidean space makes clear that this Lie group is 
not only compact and connected, but it is also simply connected [17]. The Lie algebra of 
SU{2) is usually defined in terms of the Hermitian operators Tk = \ak, k = 1,2, 3. That 
is, the basis of the representation T{SU{2)) satisfies the angular momentum algebra 

[Tk,Ti] = iekimTm, (80) 

with ekim the Levi-Civita symbol. The SU{2) group is of rank 1 (there is not a pair of inde- 
pendent elements of the algebra fl80p that commute among themselves), and the analysis 
of the structure constants = iekim shows that this is actually simple. Therefore, one 
can confine the study of SU{2) to the construction of the corresponding finite-dimensional 
irreducible representations. Remark that the set {Ti,T2,T3} is indeed an irreducible two- 
dimensional representation, as this is defined on the vector space Sp{|ef), le^)}. We are 
interested in the more general situation where the algebra flHOj) induces complex represen- 
tation spaces % with Dim('H) = n > 2. 

Let A be an n-dimensional irreducible representation of T{SU{2)) in T-t, then the 
operators A{Tk) are complex n^-square matrices. We define 

Js = A(r3), J± = [A(Ti) ± 2A(T2)], (81) 

so that 

[J+,J_]=2J3, [J3,J±] = ±J±, (82) 

and 

4 = Js, A = J^. (83) 
5.1 Irreducible representation of SU{2) 

From (!82|) one realizes that J± are raising and lowering operators for the eigenvectors of 
J3. Indeed, let \ip) be such that Jz\^) = (y\f), then 

JsJ±\y^) = ia±l)J^\^) JsJl\^) = {a±r)rj^), r G N. (84) 

Remark that the non- vanishing vectors J±\(p) are linearly independent as they belong 
to different eigenvalues of J3. Let \(pext) be such that Jslipext) = j\Vext) and J+\fext) = 
0. The eigenvalue j corresponds to the highest weight of the representation since the 
extremal state \(pext) is annihilated by J_|_ (otherwise the representation would be not 
finite dimensional). In the same context we should have J!!^\(pext) = for some positive 
integer m. Let s -|- 1 be the smallest value of m for which this is true, then one can verify 
that necessarily s = 2j. In this form, according to s, the highest weight j is either a half 
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integer or an integer. The representation A is therefore 2j + 1-dimensional and this is 
characterized by the highest weight j. For the related representation space we have 

n = Sp{\ipext), J-\iPext), J-~^Wext), J-\'^ext)}, Dim(?^) = n = 2j + 1. (85) 
Now let \(pr) := JL\'^ext)-, r = 0, 1, . . . , 2j. Using ( l82l) and ( l83l) one gets 

J+\iPr) = (2J3 + J_ J+) jrVext) = 2(j + 1 - r)\ipr-l) + J- J+|^,_l). 

After r iterations we obtain 



J+\fr) = 2 



r-1 



r(j + 1 - r) + 



IV^r-l), 



SO that 



(86) 



Then, applying (183|) and after r — 1 iterations, from (!86|) we obtain the normalization 
constant 



{iPr\(Pr) = {(Pr-l\J+\'-Pr) = r{2j + 1 - r) {ipr-l\(Pr-l) 



r!(2j)! 
(2j-r)! 



(V^extlv^ext) = C*,^- (87) 



Hereafter we shall assume {^Pextl'^ext) = 1- 

The normalized states C~}j_^\ipr=j-m) , rn = —j, ■ ■ ■ ,j, integrate an orthonormal set 
of eigenvectors belonging to J3. For these, it is customary to use the following notation 



(2j)!(j-m)! 



11 ■ ■ ■ 1 J) 



where the first entrance of | j, m) refers to the highest weight j and the second one to the 
eigenvalue m of J3 that labels the specific member of the basis we are dealing with. The 
representation space ( 185!) is then rewritten as 



n = Sp{|j, J, ), |j, J - 1), . . . , |j, -J + 1), \j, -j)} = Sp{|j, m)} 
The basis operators ( IHT]) act on this last vector space as follows 

J3\j,m) = m\j,m), 



m=j- 



]9) 



J+\j,m) = A/(j -m)(j + m + l)|j,m + 1), 
J_\j,m) = a/ (j + m)(j - m + l)|j, m - 1). 



(90) 
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5.2 Hubbard operators 

In order to express the set {J3, J±} in Hubbard notation let us introduce the change 

m^mk=j + l-k, with k = 1, . . . ,2j + 1. (91) 
Then \j,m) \j,mk) = \ j,j + 1 — k), and we can define 

:= \j,mp){j,mg\ = + l-p){j,j + 1 - g|, n = 2j + 1. (92) 
In this notation the diagonal matrix J3 reads in simple form 

n 

J3 = ^mfcX^'^ n = 2j + l. (93) 
k=i 

Using the linearity of the Hubbard operators and property f|T8l) . it is straightforward to 
verify the first of equations ( l90l) : 

n n 

fc=i fc=i 

where we have used (I9T1) . In a similar form we have the irreducible representation of the 
raising and lowering operators 

n— 1 n—1 

J^ = Y,^k{2j + l-k)X^''+\ J_ = J2Vm + ^'k)X^^''\ (94) 
fc=i fc=i 

For completeness, let us express equations (l90l) in the "/c" -representation defined in (l9Tll - 
(1921). We have 

J3\j,mk) = mk\j,mk), 

J+\j,mk) = ^/ik-l){2j + 2-k)\j,mk-i), (95) 



J_|j,mfc) = ^/k(2JTT^^\j,mk+l). 

Since the Hubbard operators X^^ are real matrices, using ( l23l) in ( !93l) and ( l94l) . it is now 
easy to verify the relationships (l83l) . For simplicity, it is also convenient to introduce the 
notation 

cl = ki2j + l-k). (96) 

Remark that 

with C| defined in ( IHTl) . Hence, the expressions (IMl) can be rewritten as 

n— 1 n—1 

J+ = ^CfeX^'^+i and J_ = ^CfcX^i'^ (97) 

fc=l k=l 

The first three lowest dimensional irreducible representations of SU (2) are reported below. 
In all cases there is agreement with the matrix representation obtained in conventional 
form reported in e.g. [lO] . 
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5.2.1 Highest weight j — 1/2. 

The lowest dimension of the representation space T-L is obtained for the weight j — 1/2. 
Thus 7^1/2 = Sp{|i, i), |i, -i)}, with Dim(7{i/2) = 2 and 



7(1/2) _ 1 ^1,1 _ 1 y2,2 _ ( \ 
— 2^2 ~ I — - 

7(1/2) _ ^1,2 _ 1 \ (1/2) _ ^2,1 _ 

5.2.2 Highest weight j = 1. 

For j = 1 we have Hi = Sp{|l,l), |1,0), |1,-1)} and Dim(?^i) = 3, with 



Ji^) = Xg'i - =10 



1 

0-1 



= V2X}'^ + V2X^'^ = 



V2 





/ 



(98) 



(99) 



5.2.3 Highest weight j = 3/2. 



For j = 3/2 wehave?^3/2 = Sp{||,3\ i\ |3 i\ |3 



2/' l2' 2/' l2' ~2/' l2' ~|)}' '^it^ Dim(?^3/2) 



4 and 



(3/2) 



3 

2^4 



1 T^2,2 1 T^3,3 3 v'iA 



/§ 

1 



\ 







1 

'2 





\ 




-IJ 



'+ 



/ \/3 \ 

2 

V3 

\ / 



J 



(3/2) y 



(100) 



6 5f/(2) X SU{2) in Hubbard notation 

In this section we solve the Clebsch-Gordan problem associated to the Kronecker product 
of two irreducible representations of SU{2). That is, we arc going to find the invariant 
subspaces that integrate the complete representation space of SU (2) x SU (2) as a direct 
sum. 



34 



6.1 The product of irreducible representations 

Let Ai and A2 be two irreducible representations of T{SU{2)) with ji and j2 the corre- 
sponding highest weights. Denote 

-Hj^ = Sp{|j„m(^-=)) |m(^-^) = -js,...Js}, s = 1,2, 

the related representation spaces. Therefore 

'Hin,n) = Sp{|ji,m(^'^)) ® \j2,m^J'^) \m'~^^^ = s = 1,2} (101) 

is the nin2-dimensional representation space of A(iiJ2) = Ai ® + ® A2, with 
rii = 2ji + 1 and n2 = 2j2 + 1. Remark that A^-^^'-^^^ is not necessarily an irreducible 
representation. As the basis of operators we shall use in each space 

4''^=A,{T,), ji^'=) = [A,(ri)±zA,(r2)], s = l,2. (102) 

Following fl9m97p . in Hubbard notation we write 

Us ris—l ris — l 

fc=l k=l k=l 

with 



^''^ - Us + l-k), c'i'^ = ^k{2j, + l-k), = 2js + 1, s = 1, 2. (104) 



We now promote the latter operators to act on the entire vector space (llOip . Thus, the 
operators 

jih,h) ^ jin) ^ ^ ^ jUun) ^ jin) ^ ^ ^ jih) (^q^^ 

correspond to the basis of SU{2) x SU{2) in the representation defined by the vector 
space 'H(jij2)- is straightforward to verify the commutation rules 

and the relationships 

'jin,n)y ^ j{n,n)^ (/l^'^A^ = J^^'^'l (107) 



In Hubbard notation the matrix representation of the operators f llOSp read as 

ni n2 

JU..2) (-1"^ + m^'^) (108) 

k=l 1=1 
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and 

ni — 1 ?12 "1 "2 — 1 

jihJi) ^ 'Sr^Jiji)-^n2ik-l)+l,n2k+e_^'Sr^ ^ ^(j2)j^n2(r-l)+i,n2(r-l)+t+l (109) 

k=l 1=1 r=l t=l 

Using Theorem Ml one can simplify the above expressions to get 



p=l 



p+n2— ".2P' 



VP'P 
n\n2 



(110) 



and 



"2(^1-1) nin2-l 

j{h:h) _ ^(ii) Tk^p,p+ji2 _L ^0'2) 



-'p+n2— n2p' nin2 



(111) 



TU1J2; _ (ii) -^p,p+n2 I „(. 

~ V ^nin2 ^ S 

p=l p=l 

where p' = \^^- As an example let us show the explicit form of some of the matrices 



fillip . In all the following examples we use = 1 and = a/2. Our results can 

be compared with those obtained in conventional form reported in e.g. 



6.1.1 Weights ji = ^ and j2 



/ 


1 


1 


o\ 











1 











1 


\o 








0/ 



(112) 



6.1.2 Weights ji = \ and ja = 1 



J. 



(Xe^'^ + Xf + Xf ) + cS^^ (X,^'^ + Xf ) + cf (Xf + Xf ) 



(1) 



/ 72 
^2 







V 



1 
1 




\ 



1 



V2 
72 
/ 



(113) 
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6.1.3 Weights ji = 1, ja 



J. 



/ 1 72 








V 









\ 

^2 

1 72 















V2 

1 

/ 



(114) 



6.2 Irreducible representation 

We want to decompose the nina-dimensional space 'H{j^,j2) into a direct sum of irre- 
ducible subspaces. Using the additivity of eigenvalues in the Kronecker product (see 
Corollary TM2.2), we realize that J1+J2 is the highest weig ht of J^^' '^'^ . Given ji + j2 - A;+ 1 , 
there are k different pairs of weights m!^^'^^ and m^-^^-* such that m^-^^^+m'^-'^) = ji+j2~^ + l- 
For k = 1 there is only one way to get m*^-^^^ + m^-^^^ = ji + ja? therefore we have a single 
state ® IJ25J2) in ^(jij2) belonging to the highest weight j = ji + j2- According 

to the discussion of Section [5TT1 |ji,ji) ® |j2,j2) is the extremal state of a vector space 
'^h+h ^ ^(ii,i2) in which the matrix basis {Ti,T2,T^} is irreducible. Denote A^^+j^ such 
an irreducible representation and write {^^3''^^"'^'', J^"'^'*'"'^^} for the basis operators. Then 

%+,2 = Sp{|j, j), J^'^+^'^^Ij, j), J^'^+^'^^+^Ij, j), . . . , fy'^''\3.3)} 

with = |ji, ji) ® |j2,j2)- In this form, Dim(Hj^+jJ = 2(ji + ja) + 1 and the vector 
space fllOll) is rewritten as the direct sum T-L^j^^j^) = T-Lj^j^j.^ ®T^r, with T-Lji C H^j^^^) such 
that Uj.+j^ nnR = Sp{|0)}. 

For k = 2 there are two different forms of solving m^^'^^ + m^^^^ = ji + j2 — I. Thus, 
the pair of states ji) ® |j2, j2 — 1) and \ji,ji — 1) ® |j2, ^2) belong to the same weight 
ji+j2 — 1 in "Hqi j2)- One of them has been already included in Tij^+j^, so that this does not 
belong to Hr. The remaining vector is the extremal state of an irreducible representation 
^ii+i2-i with I j, j-l) = ® \j2, 32 - 1) for ji > J2, or \ j, = - 1) ® | j2, ^2) 

for j2 > ji, and 

^.i+.2-i = Sp{|j, J - 1), J^'^+^'^^'^Ij, J - 1), ... , J'^''^''~'^\3,3 - 1)}. 

Hence, the vector space fllOip is rewritten as T-L^j-^j^) = 'Hj-^+j^ © "Hj^+jj-i © 'Hr', with 
Hr' n 'Hj^+j2 © ^ji+j2-i = Sp{|0)}. The procedure can be repeated at will by noticing 
that each irreducible representation is 2j + 1-dimensional with j = ji + j2,ji + J2 ~ 
1, . . . , |ji — jal ~ I5 |ji — j2\- One then arrives at the expression 

"^(jl J2) = '^jl+j2 ® "^^1+^2-1 © ■ ■ • © "^1^1-^21-1 ® "^1^1-^21 = '^(il,i2) 
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The k-th term in the above direct sum has the dimension 

4 = 2(ji + j2 + 1 - A;) + 1 = 2(ji + ja) + 3 - 2fc = m + + 1 - 2k. (115) 

The following property shows that the dimension dk of the representation space Tik is 
reduced in two units as the value of j = ji + j2 + 1 — k increases in one unit. This will be 
useful in the sequel 

dk = dk-i - 2, k>2. 



(116) 



Then, as expected, for the dimension of the entire space we have 



Dim('H(j^ 



( 2j2 + l 

XI ^fc = (2J2 + l)(2ji + 1) = nani, ji > j2 

k=l 
2ji + l 

X 4 = (2ji + l)(2j2 + 1) = nirii, > ji 

k=l 



Let jo = min{ji,j2}, then = 2jo + 1 = min{ni,n2}, and 

» . ' 

no terms 



(117) 



We now look for the operators leaving invariant the subspaces Hj in f lll7p . with 
ji +32 > j > |ji — 32\- That is, we want to construct the solution to the Clebsch-Gordan 
problem A^-'^'-'^^ — t- A'--'^'-^^^ with A^-^^'-^^^ the direct sum operator 



A {ii,i2) _ A . . ffi A ■ ■ 1 



A 



ji+i 2+i-fc ' 

fc-th term 



^Iii-i2|-i © ^Iii-j2|- 



(118) 



The first k terms in the direct sum f lllSp integrate a square matrix of order 



Zk 



X 4 = Kdk + A; - 1), A; = 1, 2, . . . , 2j2 + 1- 



(119) 



Therefore, Aj^+jj+i.^ in f lllSp is a square matrix of order 4? the first element of which 
is at the {zk-i + 1, Zk-i + 1) entrance of A'--^^'-^^-' with 

Zk-i = {k-l){dk + k), ^o:=0. (120) 
In this form, the basis elements Ja^'''^\ a = 3, ±, should read 

jUuh) ^ jin+h) ^ . . . . . . ^ J(l^'^-^'2l), J(0) := 0. (121) 

fc-th term 
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For a = 3, the Hubbard representation of the k-th term in f ll2ip is given by 

dk 



T(ii+i2+i-fc) 

•^3 



l+P, Zfc_l+P 

nin2 ' 



Zq = 0, 



where we have used fllU3p . Hence, the entire operator reads 
fc=i p=i 

In a similar form we get 

no 

JO...) = J2Y1 K\-nl-''^''-''''^\ <t'^ = v/p[2(ji+j2-A:) + 3-p]. (123) 

fc=i p=i 

Let us give a pair of examples. The following results can be compared with those we have 
reported in Sections 16.1.11 and 16.1.31 



6.2.1 Weights ji = = 1/2 



^1,1 1^3,3 
^4 — ^4 



/ 1 








o\ 


_ ( J? 
















03X1 








-1 





V 0lx3 





Vo 















J. 



J. 



(I---I) 

^2 2\> 



p=l 



/ V2 
^2 




V 

6.2.2 Weights ji = 1 and ja = | 








0/ 





03X1 


0lx3 






Ed- 



4^' 


04x2 


02x4 


4' 



p=l 



-'3 







/ - 

2 














\ 









1 
2 
















+p,4+p _ 










1 
2 






3 
2 


























1 
2 









V 
















© J3 


5l' - j'j 

— -'3 


+1) 

a 


) J3 
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J 



3 1 



E 

p=i 



X'^'^' + c. 



(1.2) 1^5,6 
2,1 ^6 



/ 73 

2 

V3 







V 



(f) 



J2x3 



\ 









J3x2 



J 



(1+1) 



J 



( 1-^ ) 



1 
0/ 



7 Clebsch-Gordan decomposition of SU{2) x 5^7(2) 



Let us consider the transformation matrix 



5 - ^ 

fc,g=l 



k,q 
nin2 



(124) 



From (I2H]) and ffTTOl) we obtain 

711712 



(125) 



p,q=l 



In similar form, equation (11221) leads to 

711712 no dk 



SJ. 



i,2fc_i+r 



t=l k=l r=l 

The (p, g)-entrance of this last matrix is then given by 



S.l 



'{31,32) 



EE 



fc=l r=l 



^ -|- I Ci2j._i+r,(j, 



(126) 



where we have used fl33l) . Given g, the double sum in (11261) is reduced to a single term 
for the labels fco G {1, 2, . . . , no} and G {1, 2, . . . , dk^} such that 



2:^0-1 + ro = q. 



Then, the equality /^^^'^^^S = SJ^^^'^^' holds whenever that 



'{31,32) 



{31) I {32) _ {31) _ ™{i2) \ c _ n 



(127) 



(128) 
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On the other hand, the matrix elements of S are obtained from ( lllip to read 



(ii j2) 



S 



Sh) a _|_ „ 



'p+n2—n2p"~'P 



(129) 



Now, departing from f ll23p one arrives at the matrix elements 



SJ 



'(iij2) 



no 

k=l r=l 



Given the double sum in the latter equation reduces to a single term for fc^, G {1,2,..., no} 
and G {1, 2, . . . , d^^ — 1} such that + + 1 = Let us take i = q + 1, then this 

last condition is reduced to fll27p with = ko and r* = tq. Therefore 



5J 



Using this last result, and 
fulfilled if 



^P,.C'n' roG{l,2,...,4o-l}- 



(130) 



with 



g + 1, one realizes that J^^''''^''S = SJ. 



r(ii.i2) 



IS 



'-'1' ^'S'r, 



(i2) 



(iij2) 



roG{l,2,...,4o-l}- (131) 



The similar procedure shows that the elements of J^^'''^'^ S = SJ^^'''^^ are conditioned to 



C(p_^„2)"^P-"2,g + Cpl\+n2-n2{p-l)'^P-'^''! ~ "^P-^+l '^k^ro^ G {1, 2, . . . , dfco !}• (132) 

Now, we use Lemma A3 of the appendix and make the change 



p = ano + l3, a = 0, 1, . . . ,ni + n2 - no - 1, /3 = 1, 2, . . . , no, (133) 
to rewrite equation 



as 



rr^ii + mf^ ' m^'^ ' "^If ^ ) = {ko + ro - a - /3 - 1)S, 



Ah) 



an2+l3, 2fcQ_i+ro 



0, 



where we have used fll04p and (11270 . It is convenient to introduce the shortcut notation 

Sp^q = Sano+i3,Zkg_i+ro = S{a,(3; kQ,ro) = S^''J[° , (134) 

so that the latter expression reads in a simpler form 

{ko + ro-a-(3-l)S'j'f = 0. (135) 

Remark that given p and q in (11340 . equation (I135p holds for all the values of ko and ro 
fulfilling (11271) . as well as all the values of a and (3 solving (I133p . Thus, for another set of 
labels ko,ro, and a, (3, fulfilling respectively (I127P and (11330 . we get 



(fco + Fo-5-/3-l)5!«~' 

ct.p 



0. 
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This last property leads to a further simplification in the notation since ko and ro can be 
written without the subindex "0", whenever they satisfy equation f ll27p . In this context, 
the nontrivial matrix elements S*^'^ are now identified by using fll34p . That is, the roots 
of the equation 

k + r-a-(3-l = 0, (136) 

with k, r, and a, (3, fulfilling fll27p and f llSSp respectively, are the labels of the matrix 
elements S*^'^ = Sp^q that can be different from zero. Using the same notation, eqs. (113 ip 
and (I132p are rewritten as 

„Ul) ck,r+l I „{j2) qk,r+l _ (iij2) qk,r qk,r j_ Jh) ck,r _ (jij2) qk,r+l fTnj\ 

'-'a+1,13 ' '^/S '-'a,f3+l ~ '^k,r '^a,^-: a-l,p ^ ^ p-1 a,l3-l ~ ^k,T K^"^ ' ) 

with rG{l,...,(ifc— 1}. To recover the conventional expression for the recurrence relations 
(I137P we use (I104p and (I123p . so that the equations to solve are given by 

v/(a + l)(2ji-a)5^^+^ + V^(2i2-/3 + l)5^f^X\ = VK2j-2fc-r + 3)^^';, (138) 

and 

v/a(2ji-a+l)^^'_\0 + V(/3-l)(2j2-/3 + 2)5^:;_i = VK2j - 2fc - r + 3) (139) 
where j = ji + j2- 

7.1 Transformation into the highest dimension invariant sub- 
space. 

Let us rewrite (I138P as 



/ (Q!+l)!(2jr-i-a)! qk,r+l j /3!(2j2-/3+l)! qk,r+l _ j r\{2j~2k-r+^)\ qk,r 

Y al{2ji~a~iy. "^Q+l,/? Y (/3-l)!(2j2-/3)! "^0,13+1 ~ Y (r-l)!{2i-2fe-r+2)! "^a,/? 

It is convenient to take 5 = a + l, /3 = /3 + l, and r = r + 1 to write 



5!(2ji-Q+l)! qk,r . / (/3-l)!(2j2-/3+2)! qk,r _ / (f-l)!(2j-2A;-r+4)! qk,r 
a\{2n-a-iy. "^0,13 + Y (/9-l)!(2i2-/3)! "^aj ~ Y (r-l)!(2i-2fc-r+2)! "^^,(3 



Now, we multiply this last equation by 



a!(2ji-a-l)!(/3-l)!(2j2-/3)! 
(r-l)!(2j-2A;-r+4)! 



to get 

2ji-5+l ^ / 5!(2ji-5)!(/3-l)!(2j2-/3+l)! ^fc,r 



(2j2-/3+l)(2j-2fc-f+4) Y (f-l)!(2j-2fc-f+3)! a,/3 

, / 2j2-^+2 / a!(2ji-a)!(^-l)!(2i2-^+l)! " r-fc.r 

Y (2ji-a)(2i-2fc-f+4) Y (f-l)!(2i--2fc-F+3)! 

2j-2k-r+3 / Q!(2ji-a)!(/3-l)!(2j2-/9+l)! cfc,r 

(2ii-o)(2i2-/3+l) V (r-l)!(2i-2fc-r+3)! 
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Thereby, we can take 



S*^'^ = const. X 



(r-l)!(2j-2A;-r + 3)! 



a!(2ji -«)!(/? -l)!(2j2-/3 + l)!^ 
to arrive at the equation 



k + r = a + (3 + 1, (140) 



2ji-a+l 



+ 



2j2-/3+2 



2j-2fc-r+3 



(2j2-/3+l)(2i-2fc-f+4) ' V (2ii-a)(2i-2fc-f+4) ~ y (2ji-a){2j2~ 13+1) ' 

After multiplying by the square-root of (2j — 2k — r + 3)(2ji — Q;)(2j2 — + 1) one gets 

k + r = a + (3 + l + {k-l). 

is fulfilled whenever k = 1. Then the equation (11401) 



We realize that condition 
reduces to 



S'j^'^ = const. X 



(r - l)!(2j - r + 1)! 



r = a + /3. 



(141) 



«!(2ji-a)!(/3-l)!(2j2-/3 + l)!' 

To fix the constant in fll4ip let us complete the binomial coefficients in the radicand. 
Thus, we take const. = ^(2ji)!(2j2)!/(2j)! so that 

n 1/2 



2ji 

a 



2j2 

/3-1 



2j 
r — 1 



r = a + /3, 



(142) 



with 



m 



{n — m)\ ml 

the binomial coefficient. Finally, it is a matter of substitution to verify that (11421) is also 
a root of (11391) for k = 1. Then, the matrix elements S*^'^ transform the representation of 

Ja^''''^\ a = ±,z, from the appropriate sectors of 'H(j^j2) into the first invariant subspace 
"Hjj+jj of the Clebsch-Gordan decomposition (II 17111 ISp . Moreover, the dimension di = 
2j + l of 'Hji+j2 is the highest of the dimensions of all the invariant subspaces Hj, ji +j2 > 
j > \ji — j2\, since dk > d^^i implies di > d2 > ■ ■ ■ > d^g (see eqs. I115l and lll6p . Next, we 
are going to derive the expressions for the matrix elements 5^'^ involving the invariant 
subspaces of dimension lower than di. 



7.2 Transformation into invariant subspaces of lower dimension. 

In this section we derive the matrix elements associated to k > 2. Let us change first 
a — J- a + 1, and then (3 f3 + 1 in ( 1139^ to get the pair of equations 

V(« + l)(2ji-a)^J^ + V(/3-l)(2j2-/3 + 2)^;^,,,_, = ^ri2j - 2k - r + 3)St'lt^ 



v/a(2ji-a + l)5t\^+i + ^/3{2j2 - /3 + 1) 5^;^ = ^r{2j -2k-r + 3)S, 



k,r+l 
Q,/3+l- 
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The substitution of this system in f ll38p leads to the recurrence relation, 

v/(« + l)(/3-l)(2ji-«)(2j2-/3 + 2) S'/^,^p 

= v/«/3(2ji - « + l)(2j2 -13 + 1) 51^,^+1 (143) 



+ [/3(2j2 - /3 + 1) + (a + l)(2ji -a)- r{2j - 2A; - r + 3)] 



fe,r 



where fc + r = Q; + /3 + l. To avoid square- roots in the coefficients of f ll43p we make also 
the change 

S':, = [«!(/3)^(2ji)^(2j, - /3 + 1)^]-^/^ 5o':I+.-i^5- (144) 
Here, S'q')^_,_^_^ and i?^'^ are to be determined while the Pochhammer symbols 

(x)„ = x{x — 1) • ■ ■ (x — n + 1) and {x)n = x{x + 1) • • ■ (x + n — 1) (145) 

respectively represent the falling and rising factorials with 

(x)o = (a;)o = 1, {x)n = {x + n - 1)„, (x)„ = {x - n + l)n. (146) 

Notice that a = in fll44p leads to (3 = k + r — 1, and necessarily E^'l,_^_^_^ = 1. After 

introducing fll44p in fll43p one arrives at the recurrence relation obeyed by E^'^^ in terms 
of the a-parameter: 

(147) 

+ [{k + r-a- l)(2j2 - A:-r + a + 2) + (a + l)(2ji - a) - r{2j - 2k - r + 3)]E'^'}, 

with /3 = fc + r — a — 1. In the solving of f ll47p we take r = 1 and proceed by induction 
on a. The lowest value a = gives rise to the expression 

<ii = -{k - l)(2j2 -k + 2) = -{k- 1)t(2j2 -{k-l) + l)i. (148) 
For a = 1 one gets 

4^2 = (fc-2)(fc-l)(2j2-A; + 3)(2j2-A; + 2) = (-l)2(fc-2)2(2j2-(A;-2) + l)2. (149) 
In general, for a = i, 

E^;l, = i-lYik - i)j{2j, -{k-i) + l)e. (150) 

Remark that i = k produces E^'q = 0, henceforth i < k — 1. On the other hand, allowing 

£ = in ( ]150p we get E^'l = 1, which is consistent with the constraint indicated after 
Eq. f ll44p . Thus, making k — £ = f3 and i = a with < i < k — 1, equation f ll50p reads as 

^a'« = (-l)"(/3)a(2j2-/3 + l)a, 0<a<A;-l, l<P<k. (151) 
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The substitution of this last result in f ll44p gives the expression of S^^^ in terms of Sq 



k-i 



Qk,l 



1/2 



-11 



(2il)a 

k-l\ (2j2 -k + 2 



S, 



a 



k,l 
0,fc 

1/2 



(152) 



nk,l 



(2ji -a + l)jj_ 

where we have used fll46p . Since S is unitary, we can fix the value of Sq']. as follows 
fc-i 



1 = E 



a=0 



s, 



k,l 
0,k 



1 + 



k-1 
1 

fe-l 

(2ji - /c + 2)^ 5 
Using the addition formula Lemma A4 



{2j2-k+2)j 
(2ii)T 



+ 





1 ^ 


\ (2j2-fc+2)— 






1 (2ii-fc+2)— 



Qfc,l 



A;- 1 



(2ii - + 2)^(2ii - + 2)^^. 



(153) 



we finally get the roots 



(2ji -k + 2)— 
(2j - 2A; + 4)^ 



(2ji)a.-i 



(2j-A; + 2)fc_i' 



(154) 



As a convention, hereafter we take the positive expression in (11541) as the definition of the 
matrix elements S^'l- Therefore, (I152p becomes 



S. 



k,l 

a, 13 



(/3)^ (2j2 - /3 + l)a(2ji - a)k^ 



(2j ~k + 2)fc_i 



1/2 



(155) 



To construct the remanent matrix elements we now use the recurrence relation (11391) . 



k,r+l 

a,l3 



l a{2ji-a + 1) k,r 
r{2j-2k + 3) "-^'^ 



' (/3-l)(2j2-/3 + 2) 
r(2j-2A;-r + 3) 



For r = 1, the straightforward calculation produces 

i,;3 ^Q-1,/3 



Qk,2 



-i\a rpk,2 Q.k,l 



with 



(/3)— (2j2 - /3 + l)»^(2ji + 13- k)k^ 
1! (A; - /?)! (2j - A; + 2),_i(2j - 2A; + 2\ 



1/2 



(156) 
(157) 
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and 



^3 =(/3-l)(2j2-/3 + 2)-a(2ji-« + l) 

-1, -a, 2ji - a + 1 1 (158) 
(3-1, -2j2 + /3 - 2 _ ■ 

To derive fllSSp we have used Lemma A5 of the appendix with 3F2 is the generahzed 
hypergeometric function 



[(/3-lM2j2-/3 + 2)T]3i^2 



3^2(0,6,0; d,e]z) = 3F2 



a, b, c; z 
d, e 



+00 

s=0 



'a)s{b)s{c)s z'' 
id)sie)s s\ ■ 



(159) 



Here we are adopting the convention 3^2(0, 6, c; rf, e) = 3^2(0, 6, c; ci, e; 1). Now we take 
r = 2 to arrive at 



f -t \a Tpk,3 f\k,2 



(160) 



The expression for 0^-2/3 can be obtained from Q'^^i ^ in fll57p after the change 1 — )■ 2. 
For the first factor in fll60p we have 



= - l)(2i2 - /3 + 2)F5_i - «(2ji - a + 1)^^^, 



A:,2 



i:m)'(^) 

s=0 ^ ^ 



»s(/^ - l)2-s(2ji - a + l)^(2j2 - /3 + 2)21 



(161) 



2)2(2j2 - /3 + 2)2 ] 3i^2 



—2, —a, 2ji — a + 1 
(3-2, -2j2 + (3-3 



In general, one can apply induction on r to verify that the roots of the system (I138m39|l 
for A; > 1 are given by 

^ir=(-i)"<r'0'-./3' (162) 

with 0^'!^^ the immediate generalization of fll57p . 

(/?)— (2j2 -(3 + l)a-r{2n +(3-kh ^ 



jk—r 



r! (fc - /?)! (2j - A; + 2)fc_,(2j - 2k + 2\_ 



and i^Qfl^^ the generalization of (116 ID . see Lemma A5 of the appendix, 



<'r' = P - Of(2j2 - /3 + 2)H 3i^2 



— r, —a, 2ji — a + 1 
(3-r, -2j2 + P-r-l 



(163) 



(164) 



To close this section we give the explicit form of the matrix S for two of the cases discussed 
in the previous sections. If ji = 32 = \ one has 

/ sl^l 



S 



s'^l 
"0" " "s^'l " "0" 



V Q s[ 



1,3 



\ / 1 

2,1 



'-'0,2 

•^1,1 

/ 
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75 



75 
\ 1 



\ 

1 

' ~i~ 

V2 

/ 



while ji = 1 and j2 = ^ leads to 



S 



•^0,1 




















q1,2 
'-'0,2 








r,2,l 
'-'0,2 








qT^ 
•^1,1 








'^1,1 











cl.3 
•^0,2 








r<2,2 
'^1,2 








'-'2,1 








r3,-2 
'-'2,1 


V 








q1,4 
'-'2,2 










7.3 Addition of angular momenta 

Consider a bipartite system integrated by independent subsystems S = 82- Any of 
the observables belonging to S should be constructed in terms of the Kronecker products 
A® with A and B observables (identities included) of S\ and S2 respectively. The 
product SU{2) x SU{2) defines the symmetry for the coupled system S = Si + S2, where 
the Lie group SU(2) represents the symmetry of each of the subsystems Si and 1S2. That 
is, the component Sk is characterized by the requirement that the n^-dimensional vector 
space Sp{\e^^)}^^^^ defines an irreducible representation of SU{2) for k = 1,2 (This is 
the situation for the and ^^C nuclei which are both SU{2) nuclear spin systems, and 
^H^^C which is a SU{2) x SU{2) system ^T]). Then, the angular momenta Ji, J2, as 
well as their sum J = Ji + J2, are conserved. These quantities are respectively the 
generators of the rotations acting on the first system alone, the second one alone, and 
both systems simultaneously [U] . The diagonalization of J ■ J and = Ji^ + J2Z gives 
rise to the orthonormal basis vectors | J, M) belonging to the eigenvalues J{J + 1) and M 
respectively. The two subsystems in turn are described by a state vector which belongs 



to Sp{|ji,m(^i)) (g) |j2,m^^'^)}, with J = ji + j2 and M 
useful to know the coefficients of the change of basis 



m 



Ui) _|_ fnih)_ This is therefore 



|J,M) = Yl {ji,m^''^;j2,m^^'^\J,M)\ji,m^^^ 



(165) 



^{n).m.fe) 



7.3.1 Clebsch-Gordan coefficients 

Let {\J, M)}'Ij^_j be the eigenvectors of the operator (87) and S the unitary matrix 
defined in f ll24p . Since S is unitary we may write 

M = (J, M\J^'''''^\J, M) = (J, M\SU^i'''^'^S\J, M) 

= (ji,m(j'i); j2,m^^'V3^''^'^|ji>^^-''^;i2,m^^'^) = m^-^'i) +m(j2), 

where \ji,'m^-'^^'-i j2,'m''-'^^) = ® \j2,fn^''''^^) and the Corollary TM2.2 was used. 

Therefore, up to a global phase (fixed as 1), we must have 

\ji,m^^'^;j2,m^^'^) = S\J,M). 
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Hence 

(J',M'|ji,m(^'^);j2,m(^'^)) = {J' , M'\S\J, M) 

defines the Clebsch-Gordan coefficient required in fll65p . Thus, the entries of the (unitary) 
transformation matrix S are associated to the couphng coefficients of SU{2) (g) SU{2). 

8 Concluding remarks 

We have apphed the Hubbard operators in the study of the Kronecker product of square 
matrices. The former represent a shorthand notation for the direct product that trans- 
forms comphcated calculations involving large matrices or a large number of factors into 
simple relations of subscripts. Thus, the Hubbard representation is compact enough to 
facilitate the study of the algebra and group properties of the observables defining a 
multipartite quantum system, no matter the order or the number of the corresponding 
matrices. In particular, we have shown that the construction of permutation matrices, 
the identification of the corresponding permutation classes of equivalence and the con- 
struction of symmetrization operators is straightforward. All the basic properties of the 
Kronecker product of square matrices have been revisited in the Hubbard representation. 
In this framework the proofs of the corresponding theorems, lemmas and corollaries are 
achieved in easy form. As an immediate application we have constructed irreducible rep- 
resentations of SU{2) by giving concrete expressions for the involved matrices in Hubbard 
notation. The same has been done for the product group SU{2) x SU{2). The solution 
of the Clebsch-Gordan decomposition of SU{2) x SU (2) by using the Hubbard represen- 
tation lead to definite expressions for the Clebsch-Gordan coefficients of the additon of 
angular momenta in terms of the hypergeometric function 3F2. Some connections can be 
found with the results already reported in [12]. In this context we like to stress that our 
results are in contraposition to the common affirmation that "it is not possible to give a 
general expression for the Clebsch-Gordan coefficients" (see e.g. |13], pp 1023). There- 
fore, the expressions derived in Section [7] of the present work give an alternative to obtain 
the Clebsch-Gordan coefficients that is different from the calculating by iteration or the 
checking in numerical tables. We hope our results have shed some light on the matter. 
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A Appendix 



Some properties of the ceiling function [■] and the floor function [-J used through this 
paper are proven. For a given real number x the ceiling function 



\x~\ = minj^; G Z : x < 2;} 



(A.2) 



yields the smallest integer greater than or equal to x. On the another hand, the floor 
function 

[xj = max{z G Z : 2; < x} (A. 3) 

gives the largest previous integer to x. 

Lemma Al. For x G M and m, n G N it is fulfilled 

(i) [x] = n if and only ifx<n<x+l 

(ii) [x + ra] = [x] + n. 

(iii) [x + nj = [xj + n. 

fiv) 



" X " 




\x/n\ 


mn 




m 



'n + 1" 




n 


m 




-m. 



+ 1. 

m LmJ 
Proof. 

(i) From the definition of ceiling function it is clear that 

[x] = n if and only if n — 1 < x < n. 

Rearranging these last inequalities we get x < n < x + 1. 

(ii) By virtue of flA.4p we may write 



[x] — 1 < X < [x] . 



Then 



|"x] + — l<x + n< [x]+n, 
and ( ]A.4p implies [x + n] = [x] + n. 
(iii) Departing from 

[xJ < X < [xJ + 1, 

one gets 

[xJ + n < X + n < [xJ + n + 1, 

and [x + nj = [xJ + n. 



(A.4) 



(A.5) 
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(iv) Consider a continuous, monotonically increasing function /. We know that if f{x) is 
an integer then x is an integer [39]. Therefore 



r/(^)i = r/(W)i- 



(A.6) 



Let us take the continuous and monotonically increasing function f{x) = —. According 
to the statement above, if f{x) = k with /c G Z, then x = km is an integer. Hence, from 
equation (]A.6|) it follows 



" X ^ 




[Ml 


m 




m 



Make x —> x/n in the last equation the proof is completed, 
(v) According to (IA.4P we may write 



n + 1 



m 



n + 1 
1 < —— < 
m 



n + 1 



m 



(A.7) 



From the right inequality we get 

n 



< 



m 



'n + r 


1 


'n + 1 




< 




m 


m 


m 



Adding 1 in both sides one has 



n 
m 



+ 1 < 



n + 1 



m 



+ 1. 



The left inequality in (IA.7I1 implies [^^] < ^ + 1 + ^, this is true if and only if 

'n + r 



n 

< - + 
m 



(A.9) 



Combining (I A. 8 1) and (1A.9I) we arrive at 

"n + 1 



n 

< - + 1< 
m 



n + 1 



m 



+ 1. 



From equation flA.Sp we see that 



n 




'n + r 


- + 1 




im J 




m 



(A.IO) 



Using (iii) in the left side of this last equation one gets the result we are looking for. □ 
Lemma A2. Let F{i,j, k, I) be a function of the indices k, I. Then 



5^$^F(.,j;fc,o<r-iHf = ^^^(^'^■'^'O^^^- 



i,j k,l 



i,j k,l 
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Proof. Let us define p — n{i — 1) + k and q — n{j — 1) + I. Note that p,q — 1,2, ... , ran 
Thus 

m n 
i,3 k,l 



m m 



nj 



= E E E F{i,k,p-n{i-l),q-n{i-l))6l 

ij p=n(j— 1)+1 q=n(j — l)+l 

^Y.^{p', q', p", q") S', = Y,F ip', p', p",p") 



m ni 



= X1 XI F{i,i,p-n{i-l),p-n{i-l)) 

i=l p=n{i-l)+l 

m n m n 

= E E ^(^' ^, ^, ^) = E E ^(^' 0^^' 4- 

i=l fc=l i,j k,l 

Where we have used Lemma Al(i). □ 

Lemma A3. Let a, f3,ni,n2 G N be such that a = 1, . . . ni — 1; /3 = 1, . . . n2, with 
ni > n2- The ceiling function satisfies the following properties 



(i) 



J. 

n2 



1. 



(ii) p' = a + 1, where p' = 

(iii) {p — n-i)' — OL, for a > 0. 



Proof. 



(i) For any x e R we have [x] = 1 if x < 1. The proof follows by noticing that ^ < 1. 

(ii) Using (i) and Lemma Al(iii) one gets 

1 r 1 

a + 1, 



Q!n2 + /3 






= a + 


[/SI 






a H 


ri2 




^2 




n2 



V 

where /3 = p — n2p' + ^2- 

(iii) The proof is similar to the one of (ii). Yet, 



an2 + /3 — n2 



n2 



a - 1 + 



A' 

ri2 



q;-1 + 1 = q;. □ 



Lemma A4. Let n,a,b e N. We have 

{a)^i {b)n (a + 6), 



E 

i=0 



(n — i)! i! n\ 



(A.12) 



where (x)7i = x{x + 1) ■ ■ ■ {x + n — 1). 
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Proof. Prom the binomial expansion we have 

(i-^)"" = E 

Note 



oo / \ 



J=0 / \m=0 / Lm 



We make the change i = / + m to get 

"^^^ (i - m)\ m\ ^^(i-m)!m! ' 

m=0 i=m ^ ' 1=0 m=u ^ ' 

On the other hand, 



n=0 



Comparing term by term the two last equations the proof is completed. □ 
Lemma A5. The generahzed hypergeometric function 3F2 satisfies 

{dUe)r 3i^2(-r, -b, c; d, -e) = ^(-1)^ ('') (6). (c)^ {d + 1)^ (e)„ (A.13) 

.=0 

with r, 6, c, (i, e e Z. The right and left Pochhammer symbols are defined as {x)n = 
x{x + 1) ■ ■ ■ {x + n — 1) and = x{x — 1) ■ ■ ■ {x — n + 1), respectively. 

Proof. Let us take r = 1, then 

sF,{-l, -b, c; d, -e) = (c?e - be) = Q (^)^ (^)^ + (^)i- 

If now we set r = 2, 

26c 6(&-l)c(c+l) 



F2{-2,-b,c;d, -e) = 1- — + 



s=0 

The proof is completed by induction on r. □ 



de d{d+ l)e(e-l) 
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